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ABSTRACT 

I.   Modified  JACOBI's  Method  [1]  for  finding  the  eigenvalues  and 
eigenvectors  of  a  Hermitian  matrix  is  a  well-suited  algorithm  for  ILLIAC  IV. 
It  is  based  on  the  idea  of  subjecting  the  matrix  to  a  series  of  orthogonal 
transformations  that  eliminate  the  off -diagonal  elements  such  that  the  matrix 
under  consideration  becomes  diagonal.   ILLIAC  IV  with  its  parallel  structure 
provides  a  tool  for  eliminating  n  off-diagonal  elements  in  one  single  sweep, 
so  that  the  whole  process  of  making  the  matrix  diagonal  becomes  very  rapid. 

II.   Modified  EBERLEIN's  Method  for  real  matrices: 

While  Jacobi's  method  is  applied  to  Hermitian  matrices,  Eberlein's 
method  [2]  applies  a  series  of  similarity  transformations  to  a  non-symmetric 
matrix  until  it  is  practically  normal.   The  resultant  normal  matrix  is  then 
reduced  to  the  diagonal  form  [2],  obtaining  the  eigenvalues  and  eigenvectors. 
The  results,  of  course,  are  best  when  the  matrix  can  be  made  diagonal. 

This  document  presents  a  brief  theoretical  background  and  a  detailed 
description  of  both  programs,  written  in  ASK,  including  the  flow-charts. 
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I.   THEORETICAL  BACKGROUND 


A.   Modified  JACOBI's  Method: 

The  classical  Jacobi  Method  reduces  a  symmetric  matrix  to  a  diagonal 
matrix  by  a  series  of  orthogonal  transformations: 


t  t 

A  x  =  cp  A  cp    ,   cp  cp  =  I . 
r+1    r  rr        rr 


Each  transformation  cp   A     cp      eliminates   two   identical  off-diagonal   elements.      It 

r  r  r 

is,  however,  possible  [l]  to  eliminate  n  off-diagonal  elements,  n  being  the 

order  of  the  matrix  A,  by  one  orthogonal  transformation.   This  can  be  achieved 

if  the  transformation  matrices  cp   are  of  the  form 

r 


cpr  =  diag  (T^T^T^ 


'  '  Tn/2} 


assuming  that  n  is  even  and 


cos  a     sin  ol 


•si 


n  OL  cos  OL 


The  matrix  A  _  will  therefore  consist  of  2  x  2  submatrices  of  the  form 
r+1 


(A  ),      _v=  (T  A  Tt)      p,q  =  (0,1,  .  .  .  n/2  -  l) 
pqy(r+l)    '  P  pq  q  r 


For  the  diagonal  submatrix 


-1- 


<v.  ■ 


2k, 2k  2k, 2k +1 


2k,2k+-l         2k+-l,2k+l 


(r)  (r) 

cos  OC        and  sin  cr   are  chosen  such  that 


X, 


(r)     /    \ 

and   sin  <X         =  -z   I  1  -  rr- 


2  .(r)    If-    "k  , 

K.      2 11  +  ~  J         sin   ^ 


where 


with 


X     -      a<r>  a(r>  Y     -   ft2    +  3^V/2 

\  ~      a2k,2k        a2k+l,2k+-l'         \  ~   l\       \/ 


t     -  2a<r> 

\  "      a2k,2kfl' 


Since    \cc_    '  \    <  x/k,   then  cos  cc    '  will  always  he   taken  positive   and   sin  or 
will  he  of  the   same   sign  as 


(r) 

tan  2CC    ; 


2a 


(r) 

2k,  2k -Hi 


2k, 2k   '   a2kfl,    2k +1 


With  this   transformation-matrix  T.     the  new  matrix    (A.  ,  )  will  he   of  the 

k  Kk'rfi 

form 


(r+1) 

2k,  2k 


(V 


r+-l 


(r+l) 

2k +1, 2k  4-1 


-2- 


After  eliminating  n  off-diagonal  elements  of  A,  the  matrix  must  be 
prepared  for  another  transformation  by  applying  the  orthogonal  transformation: 


where 


A(r+1)  =  ^A(r,l)  +*'   W*  =  I. 


V   = 


1 

0 

o      ...      0 

0 

• 

0 

0 
0 

I 

(n-2)   x    (n-2) 
identity  matrix 

0 

1 

o      ...      o 

This  permutation  shifts  the  second  row  and  second  column  into  the 
place  of  the  last  row  and  last  column,  respectively.   In  this  way,  new 
elements  are  brought  into  the  off -diagonal  positions,  and  A(r+1)  is  ready 
for  the  transformation. 


A(r+2)  ~  ^r+l  A(r+1)  Vl' 


In  order  to  subject  all  off-diagonal  elements  to  this  orthogonal 
transformation,  the  matrix  A  is  exposed  to  a  further  transformation  r 

it*  =  I, 
after  (n-2)  orthogonal  transformations  i|r  have  been  performed.  T   is  given  by 


r  = 


i2   o 


•3- 


with  I,  an  (n-m)  x  (n-m)  identity  matrix  and  I„  an  m  x  m  identity  matrix, 
where  m  is  determined  by 


n 
m  =  index  i  of  max  E   la. .  L 

The  convergence  of  A  toward  a  diagonal  matrix  WAW  is  the  fastest,  since 
TAT  rearranges  the  matrix  such  that  the  largest  off-diagonal  element  is 
eliminated  first. 

The  matrix  is  sufficiently  made  diagonal  if  the  ratio 


T]  =  E/D 
is  less  than  an  arbitrarily  small  number 


|;   |  -  10"8  (E0/Do) 


where  E  is  the  sum  of  the  squares  of  the  off-diagonal  elements  and  D  is  the 
sum  of  the  squares  of  the  diagonal  elements. 

En  and  Dn  are  calculated  from  the  original  matrix.   The  above 
value  of  |  has  proven  to  be  sufficient.   The  almost  diagonal  matrix  A  will 
be  of  the  form 


A  =  WAWt 
m 


where 


t   ,  >t 

W  =  (cp  .  . . .  Rp    ...  t9Q  Wi  ) 
N  m-1       n-1       2  TT1' 


is  the  matrix  whose  columns  are  the  eigenvectors.   Since  the  number  of  trans- 
formations is  finite,  the  resultant  matrix  A  has  the  form: 

'  m 


-k- 


A     = 
m 


\ 


To  have  bounds   for  these  eigenvalues  one   is   referred  to   Gershgorin's   theorem 
which  states   that  if  X  is   an  eigenvalue   of  an  arbitrary  n-rowed  matrix 
A  =    (a      ),      then  for   some   k,    (l  <  k  <  n), 


\k 


-   X\    <     a. 


k,l 


1  4t,k-l'         '    k,kfl'  '^n1 


For  each  k  =  1;  2,    3j  •••>    n  this  inequality  determines  a  closed  circular 
disk,,  whose  center  is  the  eigenvalue  K     and  whose  radius  is  given  by  the  sum 
of  the  absolute  values  of  the  elements  in  row  k  excluding  a,,  =  K. 

It  may  be  added  that  the  eigenvalue  problem  for  a  complex  Hermitian 
matrix  may  be  reduced  to  that  of  real  symmetric  matrices. 

Let  the  n  x  n  complex  Hermitian  matrix  A  be  denoted  by 

A  =  B  +■  i  C 

where  B  is  real  symmetric  (B  =  B  ),    and  C  is  skew- symmetric  (C  =  -C  ).   Then 
the  2n  x  2n  real  symmetric  matrix  A  given  by 


A1. 


E 


B 


has  the  eigenvalues  \n  ,  \  :   \_,  \_:  ...;  X  ,   X     and  to  each  X.   there 

1'  1'   2'     2'  '       rr  n  j 

correspond  two  orthogonal  eigenvectors 
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u  . 

3 

V  . 


and 


V  . 

J 

-u  . 

J 

If  X.   is  the  eigenvalue  of  A  then  (u .  f  iv.)  is  the  corresponding  eigenvector. 


B.   Modified  EBEKLEIN's  Method: 

This  method  can  be  stated  briefly  as  follows  [2]:   "A  matrix  A, 

which  can  be  made  diagonal,  is  normalized  by  subjecting  it  to  a  sequence  of 

similarity  transformations  A  ,  =  U  "  A  U  ,  such  that  A  ,  is  arbitrarily 
J  e-t-leee7  efl  J 

close  to  being  normal,  i.e.,  the  matrix  C=AA-AAis  arbitrarily 
to       '  '  e    e  e    e  e 

small.   Once  the  matrix  is  normal,  it  can  be  subjected  to  algorithms  like 
Jacobi's  method  to  reduce  the  matrix  to  a  diagonal  form  and,  thus,  obtain  the 
eigenvalues  and  eigenvectors  of  A. " 

The  transformation  matrices  U  are  given  by  U  =  M  P  Q  ,  where 

e  e    e  e  e 


(1) 


M  is  a  permutation  matrix  determined  as  follows: 


Let  A"  =  M1  A  M  and  C "  =  A"  A'ft  -  Aut   A",  then  M  is  chosen  such  that  each 

e     e         „  e 

2x2  diagonal  submatrix  C,  ,  has  an  element  c„,  ,  „,  of  at  least  average 

value  of  all  the  off -diagonal  elements  of  C".   For  example,  in  order  to  brin£ 

the  off-diagonal  element  cnrir,  (u  <v),  of  maximum  absolute  value,  in  the 

e.)(e.  -  e.)t. 


uv 


position  (1,2),  M  is  given  by  I,   I_  ,  where  I. .  =1  -  (e. 
*  v  '    ''      e    &      J      lu  2v'        ij         l 

Essentially  I. .  •  A  •  I. .  has  the  i-th  and  j-th  rows  and  columns  of  A 
exchanged. 


(2) 


Pe  =  diag  (T^e), 


.(e) 


'  Tn/2 


with 


M  _ 


cos  yk     sin  yR 


-sin  yk     cos  yR 


(e) 


■6- 


If  y  is  determined  by 

K. 


tan  2yk  = 


'2k-l 


,2k-l  '  C2k,2k\ 
2c2k-l,2k     /  }' 


whe 


t    .t 


re  c .  .  are  the  elements  of  the  matrix  C  =  AA  -  A  A,  and  if  cos  2v  i 
l.l  ^  k 


of  the  same  sign  as  c 


en  as  c01  ,  „,  ,   then  (c~.  n  _n  )  .    .  attains  its  maximum  value, 
D      2k-l,2k       v  2k-l,2k' (e+l) 


(3) 


Q  =  diat 
e 


q(e)   q(e)        (e) 

1  '2      '    '  '  '  ' 


3(e)l 
V2j 


with 


o(e)  ..  e(e) 
bl   "  b2 


s(e)  . 

"  n/2  " 


cos  hx    sin  hx 
e        e 


sin  hx    cos  hx 
e        e 


x  is  derived  from: 
e 


where 


with 


and 


tanh  kxQ   =  -2  •  K2(A^)/K1(A^), 


A'  =  (M  P  )   A   (M  P  ), 
e     e  e'        e  x  e  e' 


Kn  =     Z  Dn   E,   and  K_  =  Z  (D2   +  E.2  ) 
2       km  km      1       'km    km' 
k,m  k,m 


Dkm  ~  ^a2k-l,2m-l  "  a2k,2m^ 


Ekm   (a2k-l,2m  "  a2k,2m-l'" 
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It  can  be  proved  [2]  that  A  approaches  a  normal  matrix  as  e  ■+  <». 

Considering  only  real  matrices  then  for  the  practically  normal 
matrix  A,  any  diagonal  submatrix 


A 


pq 


a      a 
PP     pq 


a       a 

L  qp    q.q.j 


is  also  normal,  where  either 


a)   a   =  a       or 

pq   IP 


b)   a   =  -a       and   a   =  a 

'   pq    qp        pp   qq 


A  generalized  Jacob!  method  is  then  used  to  reduce  A  to  the  diagonal  form 
hence  obtaining  both  the  eigenvalues  and  eigenvectors  [2]. 


II.   IMPLEMENTATION 


A.   Jacobi :  Finding  Eigenvalues  of  a  Hermitian  Matrix 

1.   Storage  Scheme 

To  demonstrate  the  storage  scheme  let  us  look  at  an  example.   Let 
ILLIAC  IV;  for  the  sake  of  demonstration,  he  a  6  (six)  EE  machine. 

Assuming  all  preliminary  tests  (see  section  2.b.l-3,  P«  1*0  are 
executed,  then  BASE,  the  N  x  N  matrix  for  which  the  eigenvalues  are  sought 
will  be  of  the  form: 

BASE:   in  EE  memory 

HE:    0     1     2     3     k  5 

loo   aoi\  a02   a03  aok      a05 

a01  all/  ai2  ai3   ^  ai5 

a02  a12  /a22  a23\  a2U  a25 

a03  al3  \a23  a33/  a3^  a35 

a04  al4   a24  a3^  f^hh  %5 

a05  a15   a25  a35  W5  a55 

To  calculate  cos  OL}    sin  QL     (see  section  2.b.4,  p.  Ik)   the  matrix  is 
partitioned  into  2x2  submatrices;  i.e.,  the  procedure  ANGLE  considers  only 
the  elements 


a2k,2k>  a2k+l,2k4-l>  a2k,2k+l'  (k  -  0,  1,  2,  ...,  2  ■  1). 


-9- 


To  achieve  greatest  efficiency;  the  algorithm  calculates  cos  a,  , 
sin  a  in  pairs  of  PEs,  i.e., 


FE 


cos  ql   cos  a        cos  a    cos  QL   cos  a    cos  a 
-sin  a    sin  a   -sin  ap   sin  QL   -sin  a    sin  OL 


In  doing  so,  the  anglematrix  AJMAT  is  formed  and  unnecessary  routing  is  avoided. 
ANMAT  will  he  of  the  form: 


ANMAT:   in  IE  memory 


PE 


cos  QL    sin  QL     0        0        0        0 
-sin  QL    cos  Qn     0        0        0        0 


.n  QL    cos  QL 


cos  Qp   sin  Q?     0        0 


■sin  Oi         cos  Q      0        0 


0 

0 

0 

0 

cos  a_ 

sin  a_ 

0 

0 

0 

0 

-sin  Q.-, 

cos  a 

Since  AM4AT  is  of  this  tridiagonal  nature,  special  treatment  is  required  to 
perform  the  transformation : 

(AJMAT)t  •  BASE  ■  (ANMAT) 

BASE  •  (AB4AT)  is  multiplied  using  Knapp's  method  [3]  with  the  provision  that 
diagonals  with  all  zeroes  are  skipped.   In  doing  so,  one  row  of  (BASE  )  = 


-10- 


BASE  •  (ANMAT)  is  computed  using  only  three  multiplications  and  additions. 
Thus,  to  multiply  BASE  •  ( ANMAT),  3N  multiplications  are  required. 

In  calculating  BASE  =  (ANMAT)t  •  (BASE1)  the  algorithm  which 
performs  this  computation  takes  advantage  of  the  fact  that  BASE  is  symmetric; 
i.e.,  only  N/2  +  1  diagonals  of  (BASE1)  participate  in  the  multiplication. 
The  order  of  simultaneous  multiplications  then  becomes: 

Example:   Let  (BASE')  he  an  8  x  8  matrix. 


1 

2 

3 

h 

5 

1 

2 

3 

k 

5 

1 

.2 

3 

h 

5 

1 

2 

3 

h 

5 

5 

1 

2 

3 

U 

k 

5 

1 

2 

3 

3 

k 

5 

1 

2 

2 

3 

h 

5 

1 

Equal  numbers  represent  simultaneous  computations.   Each  diagonal  again  needs 
only  3  multiplications  so  that  BASE,  the  final  matrix,  has  all  elements 
computed  after  3((N/2)  +  l)  +  3N  multiplications  =  3  •  (3(w/2)  +  1). 

Further  explanations  are  found  in  Section  2.b.5-7>  P*  !?•   The  matrix 
(BASE')  is  now  of  the  form 


-11- 


BASE:   in  EE  memory 


EE 


a00       °       a02      a03      a0^      a05 


aii.    ai2     ai3     ai4     ai5 


a02      a12      a22       °       a24      a25 


l03     43      °      a33     a3^     a35 


a6ii     al^     a2U     a34     a44 

ai5     a15     a25     a35      °      a55 

Now  the  matrix  is  rearranged  to  bring  new  elements  into  the  (2k,2k+l), 
(2k+l,2k)  positions  (see  Section  2.b-9,  p.  21).   Then  BASE  will  be  of  the 
form: 


■12- 


(BASE')-*   in  EE  memory 
PE       0       1 


l00      a02\    a03       04      a05 


a02      a22/  "25      "12 


"03  /a33  "35     "13 


l64     tt2^    \a34     *W     °      al4 


^05     "-25     ~35  /~55     a15 


al2     ai3     ai4    \al5     aii' 

The  matrix  is  ready  for  another  elimination  and  transformation  process  as 
described  above. 


2.   Computation 

With  no  loss  in  generality,  let 

BASE:   be  the  matrix  from  which  the  eigenvalues  are 

being  calculated 
EIGV:   be  the  eigenvector -matrix 
ANMAT:  be  the  angle -matrix 
TBASE:  a  matrix  for  temporary  storage 
EPS:    an  error-matrix 

The  program  is  subdivided  into  two  parts: 

a.  the  part  that  calls  the  subroutine  (Appendix  A) 

b.  the  subroutine  itself  (Appendix  B) 

-13- 


a)   The  calling  program:   has  to  contain 

(1)  the  "DEFINE  CALL"  statement  (standard  form) 

(2)  the  matrix  containing  the  data  from  the  eigenvalues 

are  to  he  found 

(3)  the  definition  of: 

i.   the  eigenvector-matrix:   EIGV  BLK  N; 
ii.   the  angle-matrix:   ANMAT:   BLK  N; 
iii.   the  temporary  storage-matrix:   TBASE:   BLK  N; 

Note:   The  two  blocks  of  storage  for  ANMAT  and  TBASE  are  purposely  left  out- 
side of  the  subroutine ,  so  that  the  space  becomes  available  for  the  user  after 
leaving  the  subroutine. 

iv.   the  error-matrix:   EPS:   BLK  N; 

The  address  of  EPS  is  needed  if  this  routine  is  used  in  connection  with  EBEKL, 
a  routine  which  normalizes  a  matrix.   If  EIGEN  routine  is  used  alone  the 
statement  under  (iv)  is  left  out. 

v.   the  order  of  the  matrix:   DEFINE  N  =  n#; 

n  has  to  be  an  even  integer. 

Thus,  the  actual  call- statement  becomes 

CALL  EIGEN(BASE, EIGV, ANMAT, TBASE, o|EPS,N); 

The  word  "EIGEN"  is  the  name  of  the  subroutine,  and  the  user  is 
required  to  use  that  word. 

The  terms  in  parentheses  are  optional,  but  not  their  order;  that  is, 
to  make  the  call- statement  more  general,  it  must  read 

CALL  EIGEN  (<original  matrix>, 

<eigenvector-matrix>, 
•Ctemporary  storage  matrix>, 
<temporary  storage  matrix>, 
<error-matrix>|  0, 
<order  of  matrix>); 

-lk- 


The  printout  is  up  to  the  discretion  of  the  user.   If  he,  however, 
wants  the  original  matrix  BASE  printed,  the  sequence  of  instructions 
accomplishing  that  task  must  appear  "before  he  calls  the  subroutine,  since 
BASE  is  changed  during  the  computation  and  will  contain  the  eigenvalues  on 
the  main  diagonal  after  leaving  the  subroutine. 

The  subroutine  itself  contains  one  print  statement  originating  from 
the  internal  procedure  "GERSH".   The  values  printed  are  the  radii  of  the 
Gershgorin  disks,  representing  the  bounds  on  the  eigenvalues. 

b)   The  subroutine  EIGEN  and  the  Jacobi  algorithm:   The  entry  point  of  the 
subroutine  is  at  card  image  111000  and  is  named  EIGEN. 

EIGEN  makes  a  matrix  diagonal  and  returns  the  eigenvalues,  their 
corresponding  eigenvectors  and  an  upper  bound  for  those  eigenvalues. 

At  the  beginning  of  the  program  the  registers  S,  R,  X  and  D,  the 
ACARs  0  and  1  are  saved  as  well  as  a  block  of  8  local  memory  registers, 
namely  $D32-$D39.   The  user  is  advised  not  to  use  $D0-$D31,  since  they  will 
be  overwritten.   For  him,  $D32-$D63  are  available,  and  are  restored  to 
their  original  state  before  leaving  EIGEN. 

Upon  entering  the  subroutine,  $C3  will  contain  the  return-address 
which  is  saved  in  .RETUR. 

$C2  contains  the  address  of  LIST,  which  in  turn  contains  the 
addresses  of  the  parameters.   These  are  stored  as  follows: 

.ADRA  contains  address  of  BASE 

.ADRB  contains  address  of  EIGV 

.ADRC  contains  address  of  ANMAT 

.ADRD  contains  address  of  TBASE 

•ADRE  contains  address  of  EPS 

.N  contains  n  as  given  by  DEFINE  N  =  njfjh 

After  executing  a  sequence  of  instructions  which  set  up  a  series  of 
constants,  the  actual  Jacobi  algorithm  is  entered,  beginning  with  a  sequence 
of  tests: 
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1.  RWSM:   serves  to  find  the  row  index  i  for 

n-1 
max  X   |a. .;  i  =  0,  1,  2,  ...,  n-1 

and  to  store  the  result  in  .MAX. 

2.  ANYR:   The  value  for  .MAX  is  passed  into  ANYR  for  an  any-row,  any-column 

shuffle,  with  any  =  .MAX.   This  procedure  rearranges  BASE  such 
that  the  row  with  the  maximal  rowsum  is  shuffled  into  the  place  of 
the  first  row. 

The  reason  for  1.  and  2.  is  to  bring  the  largest  elements  of  BASE 
into  such  a  position  that  they  can  he  eliminated  first.  Thus,  one  achieves 
a  faster  convergence  toward  a  diagonal  matrix  BASE. 

Note:   All  calls  of  a  procedure  are  done  through  $C3  by 

SLIT(3)  =  <hame  of  procedure>; 
EXCHL(3)  $ICR; 

When  entering  a  procedure  the  return  address  in  $C3  is  always  stored  in  .SAV1. 
If  the  procedure  calls  another  procedure  the  return  address  is  stored  in  .  SAV3' 

3«   COW:   The  next  test  is  a  so-called  threshold  check.   All  elements  of 

-k 
BASE  [21-1,21]  are  compared  to  a  value  BD.   BD  =  10  "  for  k  =  1,  2, 

k,    8.   If  the  above  test  is  satisfied,  i.e.,  BASE  [21-1,21]  <  BD, 

all  other  computations  are  skipped  (see  Flowchart  1:   MAINPROGBAM) , 

and  BASE  is  rearranged  by  a  2nd  row,  2nd-column  shuffle  (for  SHUFL 

see  8.  below).   Thus,  new  elements  are  brought  into  the  (21-1, 21 ) 

positions,  and  BASE  is  retested.   BD  changes  only  after  SHUFL  has 

o 

been  executed  (n-2)  times.   As  long  as  BD  >  10   ,  BASE  is  not 
tested  for  diagonalization.   The  purpose  of  CONV  is  to  eliminate 
larger  elements  of  BASE  first,  thus  further  speeding  up  the 
convergence  of  the  algorithm. 

k.      ANGLE:   If  the  test  in  CONV  is  not  satisfied,  i.e.,  BASE  [21-1,21]  >  BD, 
then  the  procedure  ANGLE  is  entered,  and  the  rotation-angles  are 
calculated;  i.e.,  cos  OL    and  sin  <X    are  formed  from  BASE  and 
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placed  as  2  x  2  matrices  into  ANMAT.   The  program  is  straight- 
forward and  can  be  surveyed  easily,  since  many  subdivisions  and 
comments  contribute  to  its  better  readability. 

Note:   However,  the  reader  should  keep  in  mind  that  a  parallel  machine  is 
being  used,  and  that  one  test  can  be  satisfied  for  one  register  in  one  EE  but 
not  necessarily  for  another  EE. 

5.  MULTFL:   Having  found  ANMAT,  the  off -diagonal  elements  of  BASE  are 

eliminated  by  a  series  of  multiplications;  i.e.,  (ANMAT)  x  (BASE)  x 
( ANMAT) -^  is  executed.   The  same  procedure  is  used  to  find  the 
eigenvectors;  i.e.,  EIGV  =  EIGV  x  (ANMAT)1 '•  EIGV  =  I  initially. 
The  method  used  to  multiply  two  matrices  stored  in  straight  format 
is  known  as  Knapp's  Method  [3]«   Its  advantage  over  the  log- sum 
method  is  that  only  w  multiplications  and  additions  are  needed  to 
achieve  the  result,  while  in  the  log-sum  method  N^  multiplications, 
and  (k+l)  x  W   additions  are  needed,  where  2  <  N  <  2   .   As  k 
increases,  i.e.,  N  increases,  the  log-sum  method  becomes  more 
inefficient.   Because  of  the  tridiagonal  nature  of  ANMAT  further 
efficiency  in  multiplying  BASE  •  ANMAT  is  achieved  by  skipping  the 
multiplication  when  all  elements  of  ANMAT  are  zero.   The  elements 
of  each  row  of  (BASE')  are  then  found,  after  3  multiplications, 
so  that  (BASE')  =  BASE  •  ANMAT  is  found  after  3  •  N  multiplications, 

6.  TRASFOS:   finds  the  transpose  of  ANMAT  by  changing  the  sign  of  sin  Q^. 

7.  SAMUL:   This  procedure  multiplies  (ANMAT)  by  (BASE').   Taking  advantage 

of  the  symmetry  of  the  resulting  matrix  only  the  first  (N/2)  +  1 
diagonals  of  (BASE')  are  considered.  The  multiplication  process 
is  done  by  multiplying  diagonals  with  diagonals. 
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Example: 


(ANMAT)  :   in  FE  memory 


PE 


00 


01 


10 


11 


(BASE'):   in  EE  memory 


22 


23 


a 


32     "33 
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a5h 


Vp 


l55 


00 


01 


02 


03 


'ok  Do5 


10 


11 


12 


13 


\k  Dl5 


20 


21 


22 


523     D24 


25 


30 


31      32 


33 


^ 


35 


'fco     D4l     V     ^3     D^     °^5 


50 


51 


Id 


52      53 


^ 


55 


First  step:   Multiply 


aOOboo    ailbil    a22b22    a33lD33    aw"V    a55t55 


Second  step:   Multiply  after  routing  d  =  1  left    (d  =  distance) 


aoiblo  °  b2l  a23b32  °  bii3  V^  °  bo5 


Third   step:      Multiply  after  routing  d  =  1  right 

PE:  0  12  3  k  5 

0  b50  a10b0l  °  bl2  a32b23  °  b34  a5*^5 

Fourth  step:      Add  result  from  steps   1-3* 


aooboo+  allbll+  a22b22+  a33b33+  VW+  a55b55+ 

a0lbl0  al0b0l  a23b32  a32b23  a^V  Sh\^ 


and   store   in  the  main  diagonal  of  BASE. 


BASE  =    (ANMAT)1    •    (BASE') 


PE: 


aooboo+aoibio 


aiibn+aioboi 


a22b22+a23b32 


a33b33+a32b23 


akkbkh+%5b5k 


a55b55+a5^5 
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After  shifting  the  matrix  (ANMAT)  d  =  1  to  the  right  (this  is  not  actually 
done  but  think  of  it  that  way  for  better  understanding),  repeat  steps  1-k 
starting  with  element  b   down  the  diagonal  in  step  1,  b   down  the  diagonal 
in  step  2  and  b  p  down  the  diagonal  in  step  3.   Repeating  steps  1-4  ((N/2)  +  l) 
times,  all  elements  of  the  symmetric  matrix  BASE  are  known  and  the  following 
pattern  of  execution  has  developed: 


3 
2 
1 


5 

h 

5 

3 

h 

5 

2 

3 

k 

5 
h 

3 

2 
1 


5 
k 

3 
2 
1 


5 
h 

3 
2 

1 


5 
h 

3 
2 
1. 


where  equal  numbers  represent  simultaneous  operations.   The  diagonal  numbered 
5  =  ((N/2)  +  l)  forms  an  exception.   When  filling  the  rest  of  the  matrix  only 
the  elements  in  PE  (N/2)  to  PE  (N  -  l)  participate,  overwriting  the  elements 
numbered  ((n/2)  +  l)  in  PE  0  to  FE  ((n/2)  -  1). 

Since  each  diagonal  is  formed  after  3  multiplications,  all  elements 
are  found  after  3  *  ((n/2)  +  l)  multiplications.   The  total  transformation 

BASE  =  (AJMAT)1  ■  BASE  •  AKMAT 

is  executed  with  3((N/2)  +  l)  +  3N  =  3((3N/2)  +  l)  multiplications  instead  of 
2  •  (N  *2)  multiplications  using  the  conventional  ways  in  multiplying  3 
matrices.   Looking  at  a  64  x  6k   matrix  the  transformation  is  executed  after 
3((3N/2)  +  1)  =  291  multiplications  instead  of  2  •  (N  *2)  =  8192 
multiplications  using  conventional  ways. 

8.   SHUEL:   To  bring  new  elements  of  BASE  into  the  (21-1, 21 )  positions  for 

another  elimination,  BASE  enters  the  procedure  SHUEL  in  which  the 
2nd-row  and  2nd-column  are  brought  into  the  place  of  the  last  row 
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and  last  column,  respectively.   This  process  corresponds  in  theory 
to  the  transformation  cp  x  BASE/  v  x  cp  . 

Now  BASE  is  ready  for  another  transformation;  i.e.,  steps  3-8  are  re-executed. 
After  (n  -  2)  repetitions  of  these  steps,  BD  changes  from  10"  to  10  l'   ,  and 

o 

the  algorithm  is  started  from  1.   After  BD  has  reached  the  value  10"  ,  the 
matrix  BASE  is  tested  for  diagonalization.   The  convergence  factor  is  found  in 

9.  ADDIT:   Calculating  first 

E  =  Z  a. .  for  all  i  and  j, 

then 

n-1 
D  =  £  a.  ., 
i=0  1X 

one  has  S  as  S  =  E/d.   S  is  then  checked  against  KSI  =10   x  E  /D  , 
where  E  and  D  are  taken  from  the  original  matrix.   If  S  <  KSI  is 
satisfied,  the  matrix  BASE  is  sufficiently  diagonalized,  and 
BASE  [1,1]  are  the  eigenvalues.   If  S  >  KSI,  the  algorithm  is 
repeated  going  hack  to  3«  or  !• 

10.  GERSH:   Having  calculated  the  eigenvalues,  GERSH  finds  their  upper 

hounds  [h]   "by 

n-1 
radii  =  Z  a. .   i  =  0,  1,  2,  3;  • • • >   n  -  1,  i  ^  j 
j=0  1J 

according  to  Gershgorin's  theorem  for  hounds  on  eigenvalues.   The 
procedure  prints  out  the  result. 

Final  Comment:   The  routines  ROUTE,  R0TAL  and  R0TAR  are  standard  parts  of  any 
major  program.   They  adjust  routing  of  registers  and  rotating  of 
patterns  in  ACARs  to  the  left  and  right,  respectively.   They  make 
it  possible  to  handle  matrices  of  sizes  N  <  6k   and,  therefore, 
make  every  program  more  general  in  nature. 

The  flow-charts  are  purposely  made  lengthy  to  aid  the  reader  and 
help  him  gain  a  better  understanding  of  the  program. 
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Comments  within  the  program  are  made  wherever  the  author  deemed  it 
necessary.  For  the  most  part  these  comments  pertain  to  groups  of  instructions, 
Since  the  assembly  language  ASK  is  highly  mnemonic,  comments  in  abundance 
would  hamper  rather  than  facilitate  readability. 
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3-   Flowcharts 


Flowchart  1.   MAINPROGRAM:  of  Subroutine  EIGEN 


Literal  Setup 
Constants  used  in 
the  Algorithms 


Read  and  Write  the 
Original  Matrix  A 


Create  Unit-Matrix 
as  Initial  Eigenvector 
Matrix  (EIGV) 


Initializations 


© 


m       (MAI2) CZEZE 


yes 


ZEZE  =  ZEZE  +  ZEZE 


Threshold 
Check 


RWSM  (A, MAX) 


MAI1 


6 


MAX  =  row 
index  i  from 
max  Z|a. . 


ij 
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Rearrange  Matrices 
A  and  EIGV  to  Bring 
Largest  Element  of 
A  in  Position  to  be 
Eliminated  First 


CONV  Checks  if  ALL 
Super-Diagonals  are 
Less  than  ZEZE 
(Threshold  Check) 


©• 


(MAI3A) 


© 


no 


ANGLE  (A,ANMAT) 


Find  Rotation  Angle 
.  Matrix  ANMAT  from  A 


MULTPL  (A, ANMAT ; 
TRASPOS  (ANMAT) 
SAMUL  (A, ANMAT) 


Elimination  and 
Transformations 
Take  Place  Here 


MULTPL  (EIGV, ANMAT) 


0 


-2k- 


ADDIT  (A,S) 


ADDIT  Finds  the 
Convergence  Factor 
S  =  E2/D2 


no 


KSI  =  E2/D2 

from  original  matrix 


yes 


KSI  =  S 


(MAI6) 


no 


Threshold 
Check 


yes 


(MAI  *0 


& 
& 


(MAI2) 


yes 


yes 


Final  I/O 

Print  Eigenvalues 

and  Eigenvectors 
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© 


SHUFFLE  (A) 


SHUFFLE  (EIGV) 


2nd  Row-2nd  Column 
Shuffle  to  Bring  New 
Elements  into  the 
A[2I-1,2I]  Position 
for  Elimination 


(MAI3A) 


© 


Go  Back  and 
Start  a  New 
Transformation 
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Flowchart  2.   Flowcharts  of  the  Individual  Routines 


a)   RWSM:   Summing  up  the  Individual  Rows  of  A  and  Find  the  INDE  of  the  Row 
with  Maximal  Sum  (A[l,l]  does  not  participate). 


(RSU1) 


yes 


(begin  J 


Save  Return 
Address 


1  =  0 


LOAD  ROWI 


1 


Take  ABS.VAL. 
where  necessary 


Sum  up  in 
LOG  SUM 


Store  Result  in 
PEI  of  TEMP 


1  =  1  +  1 
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Find  Maximal 
Value  of  TEMP 


Compare  Register  $A 
Containing  the 
Individual  Sum  against 
Register  $S  and 
Overwrite  $A  where  $A<$S 


This  comparison  is  done 
in  the  LOOP  "RSM" .   $S 
is  shuffled  in  powers 
of  two.   Six  compari- 
sons will  fill  $S  with 
the  maximal  value . 


Compare  the  Largest 
Value  in  $S  against 
TEMP  Containing  Individual 
Sum  to  Find  Row  Index  I 


I 


Store  Index  I 
in  .MAX 


Return  to 
MAINPROGRAM 
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b)   M2i:   Stands  for  ANY-row,  ANY-column  Shuffle,  where  MY  =  .MAX.   It 

"brings  the  Row  with  the  Largest  Elements  into  a  Position,  where 
they  can  be  Eliminated  First. 


Save  Return 
Address 


I 


1  =  0 


w 

Route  ROWI  d  -  N  -  MAX 
to  bring  the  elements  of 

ROWI  into  their 
Correct  Position  within 
that  Row 

(AY) 

y 

\ 

I 

1  =  1  +  1 

BS                   y 

J  <  N 

-  1  ? 

nc 


I  =  .MAX 
J  =  0 
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(AY1) 


<>♦ 


Load  ROWI 
and  Store  it  in  J 


no 


JL 


es 


RETURN 
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ANGLE:   ANGLE  Calculates  the  Transformation-Angles  and  Creates  the 
Transformation  Matrix. 


Save  Return 
Address 


Create  Turn  On  Pattern 

and  Find  Correct 

Index  for  the  Off -Diagonal 

Elements 


Load  A[ 2x1-1, 21] 


2  x  A[2xl-l,2l] 


Distribute 
2  x  A[2xl-l,2l]  to  their 
Neighboring  Positions 
so  that  Every  PE  Contains 
that  Value  in  Groups  of  2 


Store  Result 
in  TEMPI 


0 
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0 


Load  A[I,I] 


Form  A[2I-1,2I-1]-  A[2I,2l] 


Distribute 
as  above 


Store  Result 
in  TEMP2 


no 


(AMI) 


© 


:tempi  <   o  ,      >™ (AN2I_ 


0 


32- 


0 


ANMAT    [21,21-1] 
=  -ANMAT    [21-1,21] 
=  -cos  cc     =  sin  a. 


no 


yes 


ANMAT  [21,21-1] 
=  -ANMAT  [21-1,21] 
=  cos  OL     =   sin  GL 


TEMPI  =  TEMP1/TEMP2 


O 


-0 


i, 


ANMAT    [21,21]    = 
cos   CC     = 

K 


0.5  +  o^M  tempi;     +  l 


i 


;in  ol     =s/0.5   -  0 .5/7(  TEMPI)2  +  1 


We  have  to  realize 
that  these  tests 
are  based  on  the 
parallel  structure 
of  ILLIAC  IV,  i.e. 
the  elements 
ANMAT  [21,21-1]  are 
never  the  same  for 
the  different  steps 
and  that  all  branches 
are  satisfied  for 
the  set  of  PEs  under 
consideration. 


© 
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<^TEMP1  <  0  ? 

nc 

i 

(AN3) 

sin  cc  = 
k 

-sin  oc 
k 

yes 

ANMAT  [21,21-1]  = 

-ANMAT  [21-1,21] 

=  -sin  cc, 
k 

\ 

' 

RETURN 
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d)   MULTPL :   Multiplication-Routine  A  x  B  for  Two  Matrices .   Addressing  is 

done  with  reference  to  the  location  of  BASE.   Knapp's  method  is 
applied. 


0 


© 


Save  Return 
Address 


1  =  0 


Load  Row  I 
of  A 


J  =  0 


Pick  up  Diagonal  (j) 

of  B,  Starting 

with  the  Main-Diagonal 


0 
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0 


©■ 


(MULl)   no 


Multiply 
A [I]  and  B[j] 


Form  Partial 
Sum 


Store  Result 

from  above 

Multiplication 

back   into  A 


G> 


(MUL) 


no 
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a)      SHUFL :      This   Part  brings   the   2nd-Row  to  the  Bottom  of  the  Matrix   and  the 
2nd-Column  to  the  Far  Right   and  Adjusts   the  Rest   of  the  Matrix 


ISHl 


Save  Return 
Address 


Save   2nd-Row 


1   =  1 


J  =  I  +  1 


Load  Row  J  and 
Store  it 
in  Row  I 


1  =  1  +  1 
J  =  J  +  1 


no 
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Store  Saved  2nd  Row- 
in  Row  I 


Skew  the  Matrix 


I 


Save  Second  Column 


I 


1  =  0 


I 


Route 
only 

Row  I  d  =  1  left; 
N  -  2  elements 

(SHZ) 

to  the  right  of  column  2 
participate  in  this  part 

J 

1  =  1  +  1 

no 

1   =  N-   IE 

yes 


Bring  Saved  2nd  Column 
into  Row  I 

i 

Unskew  the  Matrix 

1 

RETURN 
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)   SAMUL:   Multiplying  Two  Matrices  BASE  =  (ANMAT)1  •  ( BASE J1 


Save  Return 
Address 


I 


1  =  0 


0 


Load  NUMB 
into  Register  X 


Load  ANMAT 
Indexed  by  Register  X 
into  Register  A 


yes 

■v^-L  —  I 

no 

1 

r_ 

Route  A  d  =  I 

Right 

1 

i 

Save  Register  A  in 

TEMPI 
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Load  Register  A 
from  Register  X 


res 


Route  Register  A 
d  =  I  Right 


Store  in  Register  S 


Load  Register  A 
from  BASE, 
$S  Indexed 


Multiply  Register  A 
by  TEMPI 


Add  TEMP  to 
Form  Partial  Sum 


Store  Result  in  TEMP 


6 
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Load  $A  from  NUMB 


I 


Route  d  =  1  Right 


Store  in  $X 


Load  $A  from  ANMAT , 
$X  Indexed 


I 


Route  d  =  I  -  1  Right 


Store  in  TEMPI 


Load  $A  from  NUMB 


I 


Route  $A  d  =  I  -  1  Right 


Store  in  $S 


Load  $A  from  BASE: 
$S  Indexed 


Multiply  $A  by  TEMPI 
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Add  TEMP 
and  Store  in  TEMP 


I 


Load  $A  from  NUMB 


Route  d  =  1  Left 


Store  in  $X 


I 


Load  $A  from  ANMAT, 
$X  Indexed 


Route  $A  d  =  I  +  1  Right 


Store  in  TEMPI 


Load  NUMB  into  $A 


Route  d  -  I  +  1  Right 


I 


Store  in  $S 


I 


Load  $A  from  BASE, 
$S  Indexed 
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Multiply  $A  by  TEMPI 


I 


Add  TEMP  to  $A 
and  Store  in  TEMP 


Load  $A  from  NUMB 


Route  d  =  I  Right 
and  Store  in  $S 


I 


Store  TEMP  in  TBASE 
$S  Indexed 


I 


Load  $A  NUMB 


I 


Route  d  =  N  -  I  Right 


I 


Store  in  $S 


Store  TEMP  in  TBASE, 
$S  Indexed 


0 


no 


1=1  +  1 


_yes_ 


-H   END 
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B.   Eberlein:  Normalizing  a  Matrix 


1.   Storage  Scheme 

The  storage  scheme  is  the  same  as  described  under  Jacobi  (see 
Section  II. A. l).   The  transformations  here  are  somewhat  modified  since  we  are 
dealing  with  non-symmetric  matrices.   The  transformation  (BASE')  =  BASE  •  ANMAT 
follows  the  same  pattern  as  under  Jacobi,  however  the  transformation 
(ANMAT)t  *  (BASE1)  =  BASE  is  done  in  the  following  way: 


(ANMAT)  :   in  EE  memory 


EE 


coshi)/,         sinhty 
sinhh         cosh\[r 


coshijr  sinh\J/ 

K  K. 

sinh\(r  coshih 


cosh\[r.         sinh\J/ 
sinh\|r  cosh\|/\ 


We  notice  that  all  2x2  submatrices  are  equal  to  each  other.   Bringing 
coshij/  and  sinhnj;  up  into  CU  memory  and  then  broadcasting,  the  multiplication 
becomes : 


Step  l:  Multiply  simultaneously  cosh\h  by  row  I,  I  =  0,  2,  h,    . .  *,   W  -  2 

of  (BASE') 

Step  2:  Multiply  simultaneously  sinh\|r  by  row  I  +  1,  1  =  0,  2,  k,    . ..,  N  -  2 
of  (BASE') 

Step  3:  Add  results  from  steps  1-2  and  store  in  row  I  of  BASE 

Then  switch  coshiK  and  sinb]/  and  repeat  steps  1-3;  but  storing  the  result  in 
row  I  +  1,  now,  of  BASE.   Each  row,  thus,  is  found  after  2  multiplications, 
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(ANMAT) 


t 


achieving  a  2  •  N  multiplication  for  the  transformation  BASE 
[BASE1)-   The  total  amount  of  multiplications  for  BASE  =  (AUMAT) "    •  BASE 
UMAT  is  2N  +  3N  =  5N.   For  a  6k   x  6k   matrix  this  means  320  instead  of  the 
3192  multiplications  necessary  for  conventionally  multiplying  three  matrices 
>y  each  other. 


I.      Computation 

The  program  is  subdivided  into  two  separate  parts : 

a.  the  part  that  calls  the  subroutine 

b.  the  subroutine  itself 

The    "call"  statement  reads 

CALL  EBERL(BASE,CMAT,ANMAT,TBASE,EBEIG,EPS,N); 


AT  -  AT  •  A  (A  =  BASE) 


there 

EBERL  is  the  name  of  the  subroutine 

BASE:   the  original  matrix,  which  will  be  returned  normalized 

CMAT:    a  temporary  matrix  used  for  storing  A 

ANMAT:   a  temporary  matrix  used  for  storing  the  transformation 

matrices  P  and  Q  . 
e      e 

TBASE:   a  temporary  matrix  used  for  scratch 
EBEIG:   the  matrix  which  returns  the  product 

m    -i  -i 

e=l 

needed,  if  EBERL  is  used  in  connection  with  Jacobi, 

to  produce  correct  eigenvectors 
EPS:    the  matrix  which  returns  the  error  made  when  making  BASE 

symmetric  after  BASE  is  normalized 
N:      the  order  of  the  matrix 

lNMAT,  TBASE  and  CMAT  are  available  to  the  user  after  leaving  the  subroutine, 

The  Subroutine  EBERL  and  the  Eberlein  Algorithms.   The  entry  point 
)f  the  subroutine  is  at  card  image  155300  and  is  named  EBERL. 
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EBERL  normalizes  a  matrix  and  returns  a  symmetric  matrix  if  all 
a   =  a   or  an  asymmetric  matrix  if  some/all  a   =  -a   and  a   =  a  .   It 

pq   qp  pq    qp    pp   qq 

also  returns  the  product  of  all  transformation  matrices  and  the  error- 
matrix  found  when  making  the  practically  normal  matrix  symmetric. 

At  the  beginning  of  the  program  the  registers  S,  R,  X  and  D,  and  the 
ACARs  0  and  1  are  saved  as  well  as  two  blocks  of  8  local  memory  registers, 
namely  $D32-$Vf.   The  user  is  advised  not  to  use  $D0-$D31,  since  they  will  he 
overwritten,  unless  he  saves  the  content  of  those  registers  himself. 

Upon  entering  the  subroutine,  $C3  will  contain  the  re turn- address 
which  is  saved  in  .RETUR. 

$C2  contains  the  address  of  LIST,  which  in  turn  contains  the 
of  the  parameters.   These  are  stored  as  follows: 

.ADRA  contains  address  of  BASE 

.ADRCM  contains  address  of  CMAT 

.ADRC  contains  address  of  ANMAT 

.ADRD  contains  address  of  TBASE 

•ADRE  contains  address  of  EBEIG 

•ADRF  contains  address  of  EPS 

. N  contains  n,  the  order  of  the  matrix. 

After  these  initializations  and  after  setting  up  the  constants  pertainii 
to  the  algorithms,  the  following  procedures  are  executed  in  this  order: 

1.  MLTRPS :   This  procedure  accomplishes  the  multiplication  of  A  •  A 

(A  =  BASE)  without  actually  transposing  the  matrix  A. 

2.  TRPS:   creates  A1,  to  "ready"  A  for  the  multiplication  of  A1  ■  A.   1.  and 

2.  together  enable  us  to  calculate  CMAT  =  A  •  A*  -  A*  •  A.   Having 
calculated  CMAT  we  next  find  the  largest  off-diagonal  element  in 

3-   FWDMX:   which  returns  this  element  to  the  local  memory  with  address  . G. 
The  location  of  this  element  within  CMAT  is  also  returned  in 
.MAXR  =  rowindex  and  MAXC  =  columnindex. 

h.      AVERG :   then  tests  if  the  largest  off -diagonal  is  within  the  range  of  the 
averaged  off -diagonal  elements.   If  so  we  enter 
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5.  NULCHK:  which  determines  whether  or  not  any  of  the  c         _  c 
2k-l,2k-l    2k,  2k 

are  zero.   If  this  test  is  satisfied, 

6.  SHFT:   is  executed.   This  procedure  exchanges  the  second  row  for  row  I, 

where  I  =  row  index  for  which  sign  c   =  sign  (-c   ). 

Procedures  5.  and  6.  cause  the  difference,   con  _  ■  .  _  -  con  _.  ,  to  be 

2k-l,2k-l  2k,2k 

unequal  to  zero  and,    thus : 

2  ?   1/2 

h   =      k   r  +    (c  -    C  ) 

2k-l,2k        v   2k- 1,2k- 1         2k,2k; 

becomes  a  maximum  which  in  turn  causes  c   _  _,  =  (l/2)  h  to  be  maximal. 

2k-l,2k 

Tests  have  shown  that  including  these  procedures  avoids  oscillation  of  the 

elements  of  CMAT  when  BASE  is  close  to  normal  but  not  yet  normal  enough  to 

consider  the  results  final.   If  the  checks  under  h.    or  5-  are  not  satisfied, 

enter: 

7.  SHUT1:  which  rearranges  CMAT  according  to  MAXR  and  MAXC  found  in  3.;  i.e., 

it  brings  the  larges  off-diagonal  element  into  the  (2k-l,2k) 
position.   This  rearrangement  corresponds  to  the  transformation 

M*  •  A  •  M. 

8.  ANGL:   This  procedure  represents  the  algorithm  used  to  find  the  transfer 

matrix  ANMAT,  a  tridiagonal  matrix  consisting  of  2  x  2  submatrices 

of  the  form  P   (described  in  the  discussion  of  the  mathematical 
e 

background  of  the  Eberlein  method) .   The  actual  transformation  of 
BASE  =  (ANMAT)   '  BASE  '  AMAT  takes  place  in  the  following  sequence 
of  procedures : 

9.  MULT PL:   (BASE' )  =  BASE  •  ANMAT  is  found.   The  scheme  is  the  same  as  that 


under  Jacobi. 
0.   TRAPOS:   finds  the  transpose  of  KNMAT  by  changing  the  sign  of  the  sin 


k" 


1.  MUISA :   BASE  =  (ANMAT)   •  (BASE').   (For  a  detailed  description,  the  reader 
is  referred  to  the  discussion  of  "storage  scheme"  in  Section  II.B.l 
and  Flowchart  U.b.e.)   Having  the  intermediate  form  of  the  matrix 
BASE,  we  calculate  from  it  in: 
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12.  HYANG:   cosh   and  sinh   and  again  form  a  tridiagonal  matrix  using  the 

core  space  of  ANMAT.   HYANG  also  finds  the  convergence  factor 
IT  (A)  =  FINVAL.  A  check  to  determine  if  BASE  can  he  made 
diagonal,  i.e.,  whether  tanhU  is  not  equal  to  one,  is  included. 
If  it  is  equal  to  1,  we  leave  HYANG  immediately  and  determine 
if  BASE  has  reached  its  normal  form.   If  so,  we  leave  the  subroutine 
altogether.   If  not,  we  return  to  step  one  and  repeat  the  sequence 
of  procedures  without  undergoing  a  transformation  on  BASE  using 
ANMAT  just  found.   If  tanhU   is  equal  to  1,  we  form 
BASE  =  (ANMAT)   •  BASE  ■  ANMAT  according  to  the  procedures  described 
under  steps  9-11. 

Now  a  sequence  of  tests  is  performed  to  check  the  state  of  normali- 

-9 
zation.   First  FINVAL  <  10  '  is  tested.   If  this  test  is  not  satisfied,  we 

return  to  1.   On  the  other  hand,  if  it  is  satisfied,  we  determine  whether  or 

not  the  difference  between  the  previous  FINVAL  and  the  present  FINVAL  is  less 

-12 
than  10   .   If  not,  we  again  return  to  step  one  and  repeat  all  steps.   If  it 

is  satisfied,  the  matrix  BASE  is  sufficiently  normal  and  we  enter: 

13.  SATE:   a  procedure  for  checking  that  all  a   =  a  ;  i.e.,  we  test  sign 
pq    qp 

(a  )  since  an  accumulation  of  errors  will  never  bring  the  ideal  result. 
qP 

If  sign  (a  )  =  sign  (a   )  for  all  p  and  q, 
pq        qP 

lU.   SYM:   is  executed.   There  we  make  BASE  truly  symmetric  by  calculating 

(a   +  a..)/2,  which  is  returned  to  BASE  into  both  (i,j)  and  (j,i) 

ij    J1 
locations,  after  finding  the  error  by  computing: 

(a..)  -.  -  (a..  +  a.  .)/2  and  (a..)  ,,  -  (a..  +  a..)/2 
lj'old     ij    ijy/         ji  old     ij    ji  ' 

The  last  results  are  returned  to  the  error-matrix  EPS.   If,  however, 

sign  (a  )  =  -sign  (a   ),  i.e.,  if  a  difference  in  sign  occurs  for  some 

b    pq  qp  '     ' 

a  , 

pq 

15.   ASYM:   is  entered,  which  also  averages  the  off-diagonal  elements  by 

computing  (  a. .  +  a..  )/2.   Here,  too,  an  error-matrix  is  created  as  under  13. 
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The  subroutine  is  then  left  with  a  flag  in  $C2,  to  signal  the 

lymmetric  case  (a   =  a  )  or  the  asymmetric  case  (a   =  -a  ).   From  the 
pq    qp  pq     qp 

:ontent  of  $C2,  $C2  =0  or  $C2  =  1  respectively,  a  decision  can  he  made  to  enter 
IIGEN,  the  subroutine  for  finding  eigenvalues  and  eigenvectors  of  a  real  symmetric 
tatrix,  or  perform  the  algorithm  of  finding  complex  eigenvalues  and  complex 
ligenvectors  through  the  subroutine  CEIGEN. 
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3.      Flowcharts 


Flowchart   3.      Eberlein's   Algorithm   (MAINPROGRAM) 


Initializations 


FINVAL1  =  0 
AB  =   0 


Create  Identity  Matrix 
for  EBEIG, 
the  Total  Transformation  Matrix 


0 


C  =  CMAT  =  BASE  =  A  x 
BASET  -  BASET  x  BASE 


FINDC  (A,C,N) 


FEDMX  (C,N,MAXR,MAXC,G) 


Finds  largest  off  diagonal 
element  and  loads  it  in  G. 
MAXR  and  MAXC  are  the  row- 
and  column-indices  for  G. 
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Find  Range  between 
G  and  the  C[l,j]. 


no 


NULCHK  (C,N, INDEX) 


yes 


Exchange 
2nd  row 
with  row  I 
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0 


yes 


no 


SHUF  (a,maxr,maxc,n) 


SHUF  (C,MAXR,MAXC,N; 


SHUF  ( EBEIG, MAXR,MAXC,N) 


ANGL  (ANMAT,C,N) 


Bring  largest 
element  of  C  into 
the  (2k-l,2k) 
portion  and 
rearrange  A  and 
EBEIG  accordingly. 


ANMAT  is  the 

transformation 

matrix 


MULTPL  (A,ANMAT,N) 


:raspos  (anmat,n) 


I 


MULSA  (A,ANMAT,N) 


MULTPL  (EBEIG, ANMAT,N) 


I 


The  transformation 
A  =  ANMAT(x)a(x) 
ANMATT  takes  place 
here . 


HYANG  (ANMAT,A,N) 


Second  transforma- 
tion matrix 
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jet 


V 

MULTPL  (A,AWMAT,N) 

V 

TRASPOS  (ANMAT,N) 

V 

MULSA  (A,AKMAT,N) 

1 

f 

MULTPL  (EBEIG,ANMAT,N) 
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© 


J££_s_ 


yes 


SYM   (A,EPS,N) 


0- 


w 

AA  =  FIWVAL1 

u 

AB  =  ABS    (AA  -   AB) 

1 


$C2   =  AKTUTJl 


0 


no 


ASYM   (A,EPS,N) 


A  tag  that  can  "be 
tested  to  see  whether 
A  =  BASE  was  symmetric 
or  asymmetric. 
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Flowchart  k.      Procedures   to  Eberlein's  Method 


a)      FINMX   (C,N,MAXR,MAXC,G) 


no 


F[I]   =  max     C[I,J] 


1 


MMCOL   [I]    =  J 


1   =   1+1 


G  =  max  F[I] 


I 


MAXR  =   I 


I 


MAXC   =  MAXCOL   [I] 
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b)   FIKDC  (A,C,N) 


i)   MLTRPS 


I  =  0,K  =  0 


0 


Clear  TEMP 


J  =  0 


•»  Load  Row  I  of  BASE 


Route   NUMB  d   =  J 

_J 

' 

TEMP  =  TEMP  +  A[I]    x  A[MJMB] 

< 

f 

Route  TEMP  d  =  1 

' 

' 

J  =  J   +  1 

_%es_ 


TEMP  =  one  row 
of  storage 
in  PEM 
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ii)      TRPS    (A,B,N) 


©■ 


yes 


k  = 

1 

w 

1   =  0 

Route   NUMB  d   =  I 


I 


Load  A[NUMB] 


I 


Route  NUMB  d  =  N  -   I 


I 


B[NUMB]    =  $A 


1  =  1+1 


yes 


no 


1  =  0 


A  =  B 
B  =  C 
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c)      ANGL    (AUMAT,  C,  N) 


TEMPI  =  C[2k-l,2k] 


TEMP2  =  h  x   (C[2k-l,2k])   *  2 


V 

GERSHG  = 

=   (C[2k-l,2k-l] 

-  C[2k,2k]) 

1 

' 

TEMP2   =   SQRT    (TEMP2    +   ( GERSHG)    *   2) 


SINE  =  SQRT   (0.5   -   C[2k-l,2k]/TEMP2) 


AlMAT[2k,2k-l]    =   -AMMAT[2k-l,2k]    =   SINE 


'GERSHG  =  0  T> 121— 


no 


0 


© 


yes 


0 
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COSINE  =  SQRT   (0.5    +  C[2k-l,2k]/TEMF2) 


COSINE  =  -SQRT   (0.5   +  C[2k-l,2k]/TEME2) 


no 


0 


I 


yes 


no 


COSINE  =  0 


COSINE  =  1 
SINE  =  0 


ANMAT  [2k-l,2k-l]  =  AEMAT  [2k, 2k]  =  COSINE 


(  END  J 
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a)      HYAWG    (ANMAT,  A,  N) 


D  =  A[2J-l,2k-l]    -   A[2I,2k] 


I 


E  =  A[2I-l,2k]    -   A[2I,2k-l] 


KAP2   =  KAP2    +  D  x  E 


KAP1   =  KAP1    +  D*2+E*2 


yes 


yes 
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TANH  =  -2  KAP2/KAP1 

v 

AH  =  1   -   TAEH 

D  =  l/SQRT    (1   -    TAJffiSQ) 


E   =   SQRT    ((D+l)/2) 


COSH    :=   SQRT    ((E+l)/2) 


D  =   SQRT    ((E+l)/2   -    1) 
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AIMAT   [2k-l,2k-l]  =  AEMAT  [2k, 2k]  =  COSH 


ANMAT  [2k-l,2k]  =  ANMAT  [2k,2k-l]  =  SINH 


FINVAL  :=  0.5  x  KAP1  x  (l  -  SQRT  (l  -  TAMISQ) ) 


<D 
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e)  MULSA:   Multiplying  Two  Matrices  BASE  =  ( ANMAT)   •  (BASE')- 

If  ANMAT  was  created  in  ANGL,  the  2x2  matrices  going 
down  the  diagonal  are  different;  if  it  was  created  in 
HYANG,  they  are  equal.   A  tag  in  .SAV9  notifies  us  of 
this  fact. 


Save  Return 
Address 


1=0 


ye; 


Load  cosh\|/   and   sinh\|/ 

into    .SAV2   and    .SAV3 

respectively 


no 


0 


Find  Address  of 
the  2x2  Matrices 
of  ANMAT 
Depending  on  I 


Load  cos\|/T,    sin\J/T 
into    .SAV2   and    .SAV3 
respectively 


Multiply  Row  I  of 
BASE  by   . SAV2 
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© 


Store  Result  in  TEMP 


Multiply  Row  I  +  1 
of  BASE  by  .SAV3 


Add  Result  to  TEMP 


Store  Result  in 
Row  I  of  TBASE 


no 


Multiply  BASE 
by  .SAV3  and 
Store  in  TEMP 


J 


Multiply  Row  I  +  1 
of  BASE  by  . SAV2 


yes 


Change  the  Sign 
of  Content  of  . SAV3 
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Add  Result  to  TEMP 


Store  Result  in 
Row  I  +  1  of  TBASE 
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APPENDIX  A 
The   Calling  Program 


■  67- 


BEGIN 
FILL 


128; 


DEFINE  CALL  3NAME(BPAFAMETERS)  = 
3 IF  pSIGN(pMFIELD(pNAME))    p   THEN 
EXTERNAL  pNAME;    pFI 

3EMPTY(p PARAMETERS)   PTHEN  pELSE 
BEGIN  BLOCK 
BEGIN  USE    (63) 
LIST:      DATA  p  PARAMETERS 
END; 
CLC(2); 
SLIT (2)   LIST; 
END;   £FI 

clc(3); 

slit (3)  3 name; 

EXCHL(3)   $ICR;#; 
DEFINE     MM=U##; 

BASE:     DAr.A     3-9999970707,-3.05879213585  -k; 

LATA      -2.40792235573   -3,3-28938777238   -5,(0)60; 


DATA  -3.05879213585   -4,1.00000004Cuuw, 

DATA  -3- 9^ 56 57 ^8157  -6, -5. 9^569958816  -5,(0)60; 

DATA  -2.40792235573  -3, -3- 9^ 56 57^157  -6; 
TATA        2.00000292089,-1.72122395238    -4, (0)60; 

data      3.28938777238  -5, -5-9I+ 56996816  -5; 

DATA  -I.72122395238    -4,2.99999997304,(0)60; 


EIGV:   BLK 


\'i'p       the  eigenvector  matrix 


ANMAT:    BLK       MM;f0        THE  TRANSFORATION  MATRIX 
TEASE:    BLK       MM;%        TEMPORARY  STORAGE 
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START :      FILL; 

loll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll 'll '/l/l 'll 'll 'll '/III 'll 'll 'll 'II II II II II II II II II II II II II II II II II II II II II II II II II II 
LIT(O)  =1,3,0;%       print  out  the  original  matrix 

LIT(l)  =1,  BASE +3 > BASE; 

MAI:  DISPLAYR     $Cl,l6; 

LIT(2)        =64; 
CADD(l)      $C2; 
CROTR(l)        2k; 
CADD(l)  $C2; 

CROTL(l)        2k; 
TXEFM(O)        ,MAI; 

mil  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  imm 

CALL  EIGEN(BASE,EIGV,ANMAT, TEASE, 0, MM); 

ioii  ii  ii  ii  ii  ii  ii  ii  ii  iii/ ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  1/ ii  ii  ii  mm 

LIT(O)    =  1,3,0;%    PRINT  THE  DIAGONAL  MATRIX 

LIT(l)    =1,BASE+3,BASE, 
MAIE:    DISPLAYR  $Cl,l6; 

LIT(2)    =64; 

CADD(l)    $C2; 

CROTR(l)   24; 

CADD(l)    $C2; 

CROTL(l)   24; 

TXEFM(O)        ,MAIE; 
loll  II II II II II II II II II II II II II II II II II II II II II II II II IHI 

LIT(O)  =1,3,0; 

LIT(l)      =1,EIGV+3,EIGV;%    PRINT  THE  EIGENVECTOR  MATRIX 
MAIE1:    DISPLAYR   $Cl,l6; 

LIT(2)     =64; 

CADD(l)    $C2; 
CROTR(l)  2k; 

CADD(l)    $C2; 

CROTL(l)    2k; 

TXEFM(O)    ,MAIE1; 

ioii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  li  ii  m 

END      START. 
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APPENDIX  B 
The   Subroutine  EIGEN,   Jacobi's  Method 
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****** 

DEFI 


****** 

SDEFIN 

X 

% 

X 

X 

X 

X    . 

X 

X 

X 

X 

X 

X 

% 

X 
X 
X 
X 

Xf«l 
Xfti 

?**! 

%ti> 

Its; 
Xif  i 
Xfti 

Xffi 

?**; 
Xif; 


** 


** 


.ONE  I 
tNi 

,NMOi 

•  EnBi 

,spec» 

.ROUT  I 

.MAXi 

.A0RF.S 

.ADRES 

.ADRES 

•  KSI  l 
•ANTuN 

,antun 

•BOUNU 


BEGIN 
FILL 
tttnttitt 

NE  CALL  * 

UF  4S 

EXTERN 

SlF  &E 

BEGIN 

REGI 

LI 

FNOl 

CLC( 

SLIT 

ENn>  & 

CLc(3) 

SLITC3 

E*CHL( 

I       wRjREM 
CLCC1  ) 

cAonc l 
cAnod 

CSmL(  1 
CAnoc 1 
CCr(1  ) 
CLc(O) 
CAoOCO 
CShLCO 
CCq( 0) 
DlS^LA 
LlT(2) 

cAnrx  l 

CRnTRC 

CAnn( l 

CRqTLC 
TXpFM( 

since  in 
rase  , 

ElGVj 
ANMATt 
THASF  t 
ALL  THES 
AS  PARAM 
staustautt 

EaM 
EQij 
EUlJ 
FU|J 
E«U 
E«n 
EUli 
EUll 
E«n 
E«U 
E«U 
E(in 
EQD 
EQil 
EUU 


*#**     CALL 
NAME( 1  oARA 
I  G  N  (  K  M  F  I  E  L 

AL      ^name; 

MPTY(*PAR« 
8L0C  < 

N    USE     (63) 
ST!        nATA 


•    DEFINITION    **#***«#******<******«*##*<*#** 

METERSJa 
UUfJAME))     *TriEN 

*FI 
METERS)    UHEN    &ELSE 


&PARAMETERS 


2); 

(2)    list; 

FI 

I 

)       *  ma*e; 

3  )       «  ]  cR  »  **  '• 

t    END    nEFlNE    #********##*#*8*«*#*****### 


YR 

) 

1  ) 
) 

1  ) 
0) 

ifti 

DIRE 

THE 

ARE 

THE 

THE 

A     TE 

E     MA 

ETER 

8  *  «« 


i 

1  ! 

2: 

l 

li 

l 


SC3  J 

Sri  3  J 

?a; 

*c3; 

IS! 

Sn3  J 
2«  5 
is; 

sci »-l 
3*u ; 

Rr2J 

?/» » 
sr?; 
2a; 
»-9;«* 

r  r    atk 

MATRIX 
R  FJ  I  'N  G 

Etgenv 

A  'j  r.  L  E 
■1P0RA.V 
TRICES 
S    rRUH 

«««**» 

«,0n  »  i 
SOl  J* 
%n2i  * 

Sn3;  i, 

%n5  J  4 
So*  s  * 

*nsJ  * 
inin 
snioj 
$ni  11  * 

S0l3> 
$0l4;< 


6) 


I 

ESSTN1    IS    U^EO    LET 

FROM    rtHlCH    THE    EIGENVALUES 
FOUND. 
ECFIRMATRIX 
MAT"I  < 

Y     STORAGE    MATRIX 
NEED    TO    t)E    BY    ADRESS    PASSED 
THF    OUTSIDE 
89ttttt8**mtn*Mttititia**i*tttttiutit$»n 

FIXEU    PnlNT    Ze»t 

FlXEU    POINT    ONE. 

ORDER    OF    MATRIX 

N-l  5 

lUOOOOOOOOOOOOOnOOOOOOt 8»ENA8i ING  ONE  PE 

ENABLING  PATTERN  FOR  THE  FIRT  n  PE  S. 

64-n»  CONSTANT  USE  D  IN  END  AROuNn  ROUTING 

row  index  for  max.  val.  found  in  Rwsm. 

AdoR^SS  SAVED  HERE 


CONVERGENCE  FACTUR 

TURN-ON  PATTERN  For  THE  ANGLE  ROUTINE 

BOUND  OF  CONVERGENCE  FOR  InTERmED*  CASE. 


0000010  i 
00000?01 
0000030'' 
0000000  ") 
00000501 
0000060T 
0000070^ 
0000080  i 
00000901 
0000100  1 
0000110" 
0000120^ 
0000130" 
0000140' 
0000150- 
0000160  ■: 
0000170"> 
0000150  ' 
0000190  •' 
0000200" 
0000210^ 
0000220" 

0000230" 
0000240^ 
0000250^ 
0000260  i 
0000270  1 
0000280<"> 
0000290^ 
0000300  * 
00003 10" 
0000320" 
0000331- 
0000340- 
0000350  ■< 
0000360" 
0000370  ^ 
00003811 
0000390'i 
0000400'-. 
0000410^ 
0000420" 
000043O1 
00004401 
00004501 
0000460) 
00004701 
0000480"* 
0000490-* 
00005001 
00005100 
00005201 
00005300 
00005401 
00005501 
0000560  t 
00005700 
00005B00 
00005900 
000  06  0  0  0 
00006101 


■71- 


.BDi 

.bdi  i 

,B02i 

,B03l 

.INOE* 

.St 

.CONVt 

.SAV1« 

,SAV?I 
,SAV3» 

.SAV4i 

,SAV8I 

.SAV5« 
.SAV9S 
•  RETitH 

.AoRA  I 
. ADRHI 
.AORC  J 
. ApRD1 

,apre» 

.SAV6I 

TEMPI 

TEMP1« 

TEMP2I 
MAnTi 
GERSHu 
NUMBl 


MESSOJ 

TACK  I 

XMESSl 

SMESS2 

SMESSJ 

XMESS" 

ADBSAV 

%       **♦ 

ROUTE* 


EQli 
E«D 
E«ll 
EUij 
E«ij 
EOll 
EQll 
EUU 
E«ij 
E«U 
E«ll 
EQii 
E«ij 
fQil 
EQll 
E«ll 

EUii 
E«li 
FUll 
E«i| 

run 

E«n 

HLk 

HLk 

BLk 

RLk 

HLk 

DATA 

OAtA 

hAtA 

DATi 

RATA 

OAtA 

WDs 

DATA 

OAtA 

DATA 

OAtA 

*0s 


*nl5i* 
sni6; t 
sol  f  I 

$ni»» 
Sni*» 
so?u  l  i 
So?l  I  i 
So22J  * 

*n?5 » 
*n26  i 

*\\?7  l 

ln?v  j 
if)'30| 
Sn3i ;  t 


THRESHOLD  FACTOR  TO  BE  CHECKED  AGAINST 
THRESHOLD  FOR  SUPER-DIAGONALS. 


IN NEK  LOOP  COUNT  IN  MAIN  PROGRAM, 
CONVERGENCE-FACTOR  FOUND  IN  Anr)IT"ROUT 
CHECK  FflR  SllPFROlAGS  EoL  0. 
SAVE  REGISTER. 


********* 


*n32 

$034 

S036 
$0  3' 
1  :* 
1  : 
1  i 
1  i 
1  : 
o.  l  • 

17.1 
31  .  J 
«S»  4 
59.6 
"*#«  pf 

o : 
n»*«  PAS 
"»»«  THE 
"*«#  34S 
"#*#  Q4S 

lo: 
******** 


:  t 

> * 
:  ( 


RETUK  CONTAINS  THE  RETURN  AOR, 

TO  l.lN*  TO  THE  OUTSloE. 

AURESS  UF  ORIGINAL  MATRIX 

AORESS  UF  ElGENVtCTURMATR. 

ANGLE  MATRIX 

TEMP.  STORAGE  MATRIX 

AORESS  OF  THE  ERKOR  MATRIX 

OiNE  HOW  OF  SAVE-STORAGE  IN  PE-HEMOPY 


2.3. 
H.  19 
2.  33 

b.  47 
0.61 

I  CO 

E  AF 

TkA 

L  AF 

E  AF 


FILL  I  X 
STL(3) 
LDL(3) 
RTi 


LDL(3) 

L0EE1 

LDA 

L0l(3) 

ISL 

SETE 

SETEl 

CLrA  J 

LDL(3) 

RTL 

LOA 

LDL(3) 

COmPC( 


.*aV3;s 

.SA V2  J  < 

$a.0(3  ) ; 

.SPEC  J  * 

8r3» 

*PI 

hhihijsh 
NIJMH  I  t 
. SAV2;  * 
«r3j  * 
i .  a  n  [) .  e  ; 

F.AND.EJ 

.pout; * 

0(3)  1  t 

*R» 

.spec; 

3)5 


4.  5.  6.  7»8»  9.  1 ".  1  1  .  12.  13.  14.  15.  16, 
.  ?0»2W2.?3.?/i»?5.26.27.?8»29.3M 

.34.3  5.36.37,}8.39.4n.41./i2»4  3»44, 
»!,8»4  9.br>.51»52'53.bu.55.56»57.59, 
.4  2.63  J*     PE  NUMoERING 
NTAIUS  bOUND  ON  ElGtNVAL.  I  *  *  *  "  j 

TFR  ANY  row/col.  SHUFFLE  *##"» 
NSFOHMATl ONMATRI X  -  AMMAT  * » « " J 
TFR  ThE  TRANSFORMATION  «**"» 
TEP  2-nO  RUW/CoL.  SHUFFLE  #*#"J 

****************************************i 

ROUTE-ROUTINE  FOR  N"»64 

SAVE  RFTiiRN  AOOpEss. 

LOAU  ROUTING  0T STANCE, 
*   SA  lb  ASSUMED  TO  CONTAIN  THE  ELEMENTS 

TO  IE  ROUTED. 


THIs  PAPT  CLFANs  THE  FI-RST  0  PF«<; 
OF  RbA,  IN  CASE  FLEMENTS  HA\/E  "FEN 
ROJTED  INTO  THOSE  PE"S.  FROM  AccROSS 
I  ThE  ouundary  N. 


BRING  THOSE  ELEMENTS  WHICH  wERF  ROUTED 
ACROSS  ROUNOARY  INTO  THE  FIRST  n  pE«S 


B0^i000062OO 

0000630'! 
0000640" 

00006S0'! 
00006600 
0000-S70 
INE. 00006*0  i 
0000690" 
0000700^ 
OOOOMO" 
0000720'1 
00007300 
00007400 
0000750  i 
0000760" 
00007  700 
00007ri0o 
00007^0' 
0000800  i 
OOOOfllOo 
0000820  i 
0000«300 
000084  0"! 
00008SO  : 
0000660" 

00  00  9 70' 
0000890'' 
0000990" 
0000900  ' 
0000910-1 
0000920'"1 
0000930  . 
0000940" 
0000950" 
00009600 
0000970" 
0000990'' 
0000990') 
00010000 
0001010" 

****00010200 
0001030" 
0001040O 
0001050O 
00010600 
0001070" 
0001080" 
00010  9  0  0 
0001100"! 
0001 110" 
0001120O 
0001130  ) 

ooonaoo 
0001150"! 

0001  1^.00 
0001 1700 
0001 190" 
0001 1900 
0001200" 
0001210O 
00012200 
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tf IM$ ** 


I!   **** 
ROTARI 
%       AQJU 


1$*$*** 

ROTAL" 
%       ADJU 


ttatts* 

RWSMi 

% 

% 

% 

% 

% 


LDEE1 
ClPAJX 
COmPCO) 
LDpEl 

LDlO) 
rXCHL(3) 
********** 

f  ill  i 

ST  ACAR1  F 
STL(3) 

LDLn) 

CRnTRd  ) 

LDlC?) 
CShR(2) 
CAnH(2) 
CAmo(2) 
CExOR(  1  ) 
L0i(3) 
CRnTK  1  ) 
CEXORC 1 ) 

STL(l) 
LD|_(3) 
EXCHL(3) 

f  ill  j 

ST     AC4R1     F 
STlO) 
LDL(l) 
CRflTLC  1) 
LD|_C3) 
CShL(3) 
CAnoO) 
CAnoO) 

cExoru ) 

L0L(2) 
CShL(I) 
CEXORU ) 
STl( 1) 
LDl(3) 
EXrHL(3) 

FILL  J 


»C3 


*c3> 


CLEAR    RfiA    PAST    THE    BOUNDARY    N 


rttnnntttsttittt»tn»$t*ttittt«itt»ttttttttiuttMtt 

.sav3js    and  the  result  <iill  8e  in  sa  properly 
stcR»*     Routed. then  return. 

****************************************************** 


OH 


shift, right  for  d  lss  64 

RETURN  aorfss 
PATTERN  THAT  NFEqS  ADJUSTMENT 
ACAHQ  CONTAINS  L) 


RSUi  I 


STL(3) 

sEte 

sEtei 

clraj 

LDs 
LlT(O) 

cAon(O) 

CRpTL(O) 
LDL(1  ) 

CCr(  1) 
LDL(2i 
CAoO(2) 


AN  E\D-AROUND 
.SAV3J*   SAVE 
,  $  A  V  a ;  % 
0(0 )  :x 
.spec; 
oco) ; 

*c 1 ; 
snt>8 
*c2; 

,m; 

0(3)! 

.  S  A  v  6 ; 
,sav3; 
MCR: 

or  am  end-around  shift.left  for  d  lss  6« 

.sav3; 

.SAV-55 

o<o>> 

.SPECS 

0(0)  ; 

$rlt 

SD5S 

tC3» 

.  RU'JT  J 

0  (  2  )  i 

*C35 

.SAV3J 

.SAV3! 

MCR; 

Knz<tsitii2ttiattt«iUgsiitg*s$*ttti*g»$stB 


SUM  UP  IN  ARS.  VALUE  EvERY  RO^j  OF  THE 

MATRIX  BASE  AND  FIND  FIRST  THE  MAXIMUM 

OF  THE  INDIVIDUAL  SUMS  ANO  THEN  THF 

INDEX  OF  THE  ROW  WHERE  THE  MAvlM(jM  CAME 

FROM. 

SAVE  RETURN  ADDRESS. 


.  S  A  V  1  i  I 

E.PR.-EJ 
E. ANd.EJ 


*A  J 

=0.1»0»X    LOOP  FOR  PICKING  UP  THE  ROwS  OF  THE  MATRI 

2Ut 

.SPECJX 


0(0)  J* 
. AORAf 
SCO» 


THE  ELEMENT  A[I,I]  DOES  NOT  PARTICIPATE 
IN  THE  SUMMING. THEREFORE  WE  TURN  OFF 
PECU  BY  THIS  CONSTRUCT 


000123CO 
00012400 
0001250^ 
0001260"> 
0001270"> 
00012flOi 
O001290  ^ 
*0001300  > 
0001310  > 
0001320  "> 
0001330- 
0001340  ^ 
0001350" 
0001360  i 
0001370  - 
0001390  ■) 
0001390- 
0001400'- 
0001410"! 
0001420  > 
0001430- 
000144  0^ 
0001450  i 
0001460"' 
0001470' 
00014*10-1 
0001490  ■* 
0001500- 
0001510  J 
000152C"1 
0001530- 
000154C" 
0001550- 
0001560^ 
0001570'^ 
0001590'! 
0001590  i 
0001600  ^ 
000161 C- 
00016201 
0001630  * 
00016400 
0001650  1 
00016  6  0  ■ 
0001670^ 
0  0016  9  0--. 
0001690  i 

0001700"! 
000171CO 

0001 720- 
0001730". 
00017401 
0001750-^ 
<  0  0  0  1  7  6  C  0 
0001770-1 
0001780., 
00017900 
00018000 
00018100 
00018200 
00018300 
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LDEE1 
IDA 
CLCC  3) J 

ldeei 

LlT(l) 

LDA 

LDLd) 
LDfCl 
ADRN 
% 

LOs 

TXETMCO) 
JUMP 
STS 

sEte 

StTEl 
LlT(2) 

LOS 
LOA 
L0l(2) 

ldeei 
lda 

StTE 

setei 

LDS 

LlT(l) 
RSMI       RTL 
LOs 

ial 

StTE 
SETE1 

LOA 
SETE 

setei 

LOs 

CAnoci ) 

LIK2) 

EULXFCl ) 


t 

ANY 
t 
% 
X 


LO|  (2) 

LDEE1 

LDA 


SETC(1  ) 

compcc i ) 

CAnOC  1  ) 
LEADOC 1  ) 
LlT(O) 
CAmD( 1 ) 
STl(1  ) 

LDl<  3) 

EX^MLO) 
************** 

R«    r  ill  » 


0(2)1 

sc3» 
■  o.oj 

.SPEC 

sen 


FORM  PARTIAL  SUM.  «E  KNHW  THE 
MATRIX  IS  SYMMETRIC 


SA5 

RSill  '" 

Tfmp: 

E.TH.-EJ?   FIND  MAXIMAL  VALUE  IN  TEMP, 

E  .  A  N  1  .  E  f 

«0.0» 

sr2» 
ssj 

.S°EC I 

*C2> 
TFMPJ 
E.0H.-E5 
E.  ANQ.EJ 

sa; 
=  i  J  • 

ss.oc  i ); 
%o  t 

SSI* 

I  . A  N  0 . E  i %       CUMPARFQ  WITH  $  S  • 

E.  AN().EJ 

*SI 

e.tr.-e; 

E. ANn.EJ 
S  A  ) 

*ci » 

=  64» 

sc2»RSm; 


X   ROijTE  SS  IN  PQWFRS  OF  TWO 

ELIMINATE  SMALLFK  VALUES  OF  *A  whE 


.SPEC t 
SC2» 
TEMPJ  i 


S^J 

II 
I 

5"?i 
) 

x77»8; 

SCO' 
.MAX} i 

.SAVl 1% 

siCr» 

********** 


i   MAXIMAL  VALUE  IS  FOUND.  wHEN  TEST 
WHERE  niO  THE  LARGEST  VALUE  Come  F 
FIND  ROW  INDEX. 


at  thi?  point  every  pe  contains  th 
mavlmal  value  in  ss.  while  the  fir 
pe  s  cnntaln  the  individual  sums  i 


.MAX  CONTAINS  THE  ROW  INDEv  WHERE 

VALUt  CAME  FROM, 

RETURN  TO  ThE  MAIN  PROGRAM. 

************************************* 

ANYR  BRINGS  THE  ROw  WITH  THE  LARGE 
INTO  THE  TOP  Rn*i  ANO  PULLS  ALL  OTH 
IT  IN  A  CIRCULAR  MOTION. 


00018ft0<"> 

00018500 

00018600 

00018701 

0001880' 

00018901 

00019001 

0001910  i 

00019201 

0001930') 

00019*01 

000 1 95  On 

0001960^ 

00019700 

00019801 

00019901 

0002000^ 

0002010  ) 

00020201 

00020300 

0002040  ' 

0002050  ! 

0002060  l 

0002070> 

00020801 

0002090  ! 

0002100" 

00021101 

N        00021201 

0002130  • 

0002140" 

0002150 

00021601 

0002170  ~ 

0002180" 

0002190" 

0002200 

FAILS.   00022101 

ROM      0002220' 

0002230  ^ 

00022401 

00022501 

E        0002260: 

ST  N     0002270" 

N  $A.    0002280" 

0002290  ' 

0002300-' 

0002310'' 

0002320" 

0002330  1 

0002340" 

0002350- 

LARGEST  0002360" 

00023701 

00023801 

0002390" 

********00024  00  J 

0002410  ! 

ST  SUM   00024201 

FR  wITH  00024301 

00024400 
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STlO)  .SAV1 

L0l(2)  .Ml 

CSijB(2)  S07JX 


l  i 


STl(2)  .SAV2J 

LDlO)  .SPEC? 

LDeEI  J.C3J 

L  I T  (  0  )  =  0  .  1  • 

cRdtlco)  ?m% 

cAr)O(O)  Sn3; 

AYl       LDL<2)  .AOHE 

cAno(2)  sro; 

IDS  0(?)» 
CLC(3)f 
SLTT(3)  : 
STA 

LlT(O) 

LlT(  1  )  =0.1 

CAnn(  l )  *n3; 

c  h  n  t  l  ( 1 1  ?  'i  i 

CAnn( 1 )  *n7  s 

AY1«     L0K3)  .adR 

cAnnc3 )  %r\  i 

LDa  Of  3  ) 

l_Di_(2)  .  inR 

CAno(2)  sroJ 

STfl  0(?) 

ALIT(O)  =i; 

T^rfMll)  .AYl 

LlTfl)  =1.0 

CRnTR(l)  24J 

CAnO( 1  )  lr7 ; 

CSuR(l)  *0l * 

CRnTLU)  241 

AY2»      LDL(3)  .AnR 

CAnO(3)  ir 1  J 

LDA  0(3)! 

LDl_(2)  .  ADRl) 

CAnD(2)  *,r.ni 

STa  0(2); 

ALIT(O)  =i; 

TXeFm(I)  .AY2) 

LlT( 1 )  =1.0. 

CRDTR(l)  2'ilt- 

cAnn(i)  %^m 

cRnTL(l)  ?a; 

AY3 1      l_0l(3)  .40RE 

C AODC  3  )  SCI  ' 

tD[_(2)  .AnHD 

CAnO(2)  SC1J 

LOA  0(2)» 

STa  0(3)* 

TX^FHCl)  . aY3 J 

LDl(3)  ,<;A  VI 

EXCHL(3)  MCHJ 

j   ********************* 

SHUFL«    FlLU 


SAVE  RETURN  ADDRESS. 

N-.MAX.  MAX  IS  THE  ROw-INOEX  FOR  THE 
LARGEST  RUw-SUM.  N-MAXIS  THE  ROUTING 
DISTANCE. 


BRING  THE  ELEMFNTS  INTO  THEIR  RESPECTIVE 
POSITIONS  WITHIM  THEIR  RO-^S. 


0;% 


s; 


R0i|TE;EXCHL(3)    SICR;«    GO    TO    ROUTE    ROUTINE 
0(2>! 

.  \Y  I  t 

al)  :  * 


THf    EVENTS     ArE     In    THEIR    NEW    POTIONS 
NU*    *E    HRINC,    THE    ELEMENTS    INTO    TmEIR 
FlMAL    LOCATION     Ri    CIRCULAR    MUTION 


uJ 


es; 
i 

OS 

J 


0  1 


ES> 


o;* 


SJ 


BRING  THE  ELEMENTS  FROM  T*ASE  INTO 
THEIR  RESPECTIVE  MATRIX. 


J 

*      GO  BACK  INTO  MAIN-PROGRAM 
*********************************************** 


0002450^ 
0002460^ 
0002470' • 
0002480  • 
0002490" 
0002500^ 
0002510  > 
0002520" 
0002530'". 
0002540-1 
0002550-: 
0002560^ 
0002570" 
0002580" 
0002590" 
0002600-' 
0002610" 
0002620" 
0002630 .. 
0002640^ 
0002650'! 
0002660 ' 
0002670 
0002680 
0002690" 
0002700. 
0002710'* 
0002720" 
0002730  ~ 
0002740 
0002750  ^ 
0002760  <; 
0002770  i! 
0002780  ■• 
0002790  " 
0002800' 
0002810"!; 
0002*20  • 
0002830  i 
0002840" 
0002850" 
0002860" 
0002870 
0002880" 
0002590  ^ 
0002900  1 
0002910^ 
0002920" 
0002930' 
0002940^ 
0002950  i 
0002960^ 
0002970  i 
0002980^ 
0002990" 
00030001 
0003010^ 
0003020'' 
00030300 
♦00030400 
00030500 
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SHI 


SHHI 


SHI  I 


SH2I 


STl(3) 
LDlO) 

c*nn(3 ) 

LDl(O) 

ldfei 

loa 

sTa 

LOl(2) 
LOlCO) 
CAnDC  0) 

cAnn(?) 

E<H  XT(2) 

LDa 

STa 

CAnO(3) 

sKtp 

LOa 

STa 

LlTfO) 

cAno(O) 

cRnTL(O) 

LOL<?) 

CAoO(2) 

LOA 

STl(O) 

CLc(3); 

SLITC3) 

STa 

TXpFM(O) 
LOa 

LOLf 0) 
STLfO) 
CLc( 3)  J 
SLIT(3) 

LOS 

LOx 

LOLfl) 

LOa 

STa 

LlT(O) 

C  A  n  0  (  0  ) 

CRnTL(O) 

LDL<?) 

LlT(3) 

CExHR(2) 

STl(2) 

clc(3); 

SLIT(3) 

L0l(2) 

LDEEl 

L0l(1) 

CAnD(l) 

STL(1  ) 

LOA 

L0l(3) 

ST|_(3) 

CLr( 3)1 

SLIT(3) 

STl(O) 


i  f> 


SHUFL  BRINGS  2ND  R0«  TO  THF  BOTTOM 
AND  *N0  COLUMN  TU  THE  RIGHT  HANO  ENOt 

PICK  UP  TrF  SErONo  ROW  OF  ThE  MATRIX 

IN  USE  AND  STORE  IT  TEMPORARILY  IN  TEMP 


H : 


.SAVl J 

. 4  ORES 
Ji|l( 
.SPECS 
*co; 

Of  3  )  I 

TFMP: 

ir  3  * 
,m  t 

I'M*  ; 

m  » 

or?); 
or.  3) ; 

SOI  J 
.  sH'J 
TFMP  J 

o  c  3 >;* 

sO. I . 0 

*n3  j 
?'»; 

.  a  n  R  E  s 
trot 
0(?x: 

. SAV2! 


=rohte;exchl( 3)  *ICR» 
or?); 

.  s  H  l ;  *     S  K  F  *  "  E  N  0 
NiimhJa      ADJUST  INDEX  TO  PICK  UP  2ND  COLUMN. 

.one; 

.SAV2J 

=  RniiTL;EXCHL(3)     flCRJ 
<A» 


;  * 


E  N  0    2  N  0    R  0  V 
SkFi    MATRIX 


SH'IFFLE 

IN    ijSE    FOR    COLUMN    EXCESS 


END  ADJUST 

PICK  uP  2N0  COLUMN  ANj  STORE  It  IN  TEMP, 

END  PICK  UP 

REARRANGE  ThE  REST  OF  THE  MATRIX 


ssr* 

.  aqkes;* 
*o( 1  )  j 
Temp; * 
so. 1 »0; x 
S'i  3 ; 

2  a  I 

.  .spec; 

=140000000000000000000018; 
$03  J 
.  5  A  V  8  J 

:ROTAK JEXCHLC3 )  „ I c R J 
.SAVH; 
%C21 
.  A"iRES; 
%rni 
.SAV«; 

od  j; 

.  N-.'O  J 
,SAV2i 

Rn.iTE;ExcriL(  3)  jicr; 
.SA^b; 


000  3  0  6  0) 

00030700 

00030800 

00030900 

0003100'' 

00031 100 

0003120  i 

00031 30i 

0003140"' 

0003150  ) 

00031601 

0003170  I 

0003180" 

0003190^ 

0003200" 

0003210  i 

0003220  ) 

0003230  > 

0003240 0 

0003250" 

0003260'' 

0003270  i 

00032*0" 

0003290^ 

0003300  1 

0003310  » 

00033201 

0003330" 

0003340'! 

0003350'i 

0003360'' 

0003370" 

0003380^ 

0003390" 

0003400" 

0003410"! 

0003  4  2  0" 

0003430" 

0003441 

00034501 
0003460" 
00034701 
0003480" 
00034901 
0003500" 
0003510" 
0003520" 
0003530" 
0003540" 
0003550'^ 
0003560^ 
000  3570" 
0003580" 
0003590" 
0003600  I 
0003610.' 
00036201 
0003630" 
0003640*" 
00036501 
000366OO 
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I 


SH3I 


X 

ANG 
X 
X 
X 


lol( 

STL< 
CLc( 
SLIT 
LOL( 
LDeE 
LDL( 
LOL( 
STa 

txef 

LDL( 
LDEE 
LOA 

ldlc 

CSu» 
STL( 

CLc( 

SLiT 

LOs 

LOx 

LDA 

LOl( 

CSuB 

ST|_C 

CLC( 

SLIT 

LDL( 

STA 

LDl< 

LDl( 

cAnn 

CSuD 

sTl( 

LDa 
CLc( 
SLIT 
STa 
E«LX 
LDL( 
EXCH 
********** 


0) 

2) 

31) 

(3) 

2) 

1 

0) 

1) 

M(0) 
0) 

1 

1) 

(1) 

1) 

3)5 
(3) 


.ONE! 

.«,av*> 


'ROTAL 
.SAV 

lC.2'> 
.SAV 
.SAV 

0(1  ) 

»  SH2 

.SPE 

$ro» 

SXJ* 

tm  ; 
•  sav 


JEXCHL(3)  SIC«» 
8J 


5; 
4; 
J 

;  f> 
c;  1 


ENn  REST-ARRAN-.E 

BRING  2ND  COL.  INTO  LAST  COLt 

FIRST  j  AOJUST  INDEX 


END  iNnEX-AOJUST 


2) 

(?) 

2) 

3)J 

(3) 

1  ) 

1  ) 

2) 

(1) 

(2) 

2) 

3)1 
(3)  = 

F(2) 
3) 

L(3) 
*  **  *  * 


!kOliTE;EXCHL(3)  ilCRI 
SA  J 
$S»  * 

tfmpi 
,  m  m  0 ; 

*  01  » 

.SAV2J 


:RGllTL  JEXCHL(3)  UCRI 
.  Af^RESS 


END  COL. "SWITCH 
UNSKLw  MATRIX  IN  i)SE 


LE» 


ANl 


f  ill  ; 


ST!  (3) 
LlT<0) 

CAOO(O) 

cRnTL(O) 

CAnn(0) 

LDlM) 

CCr(1  ) 

TXfFM(O) 

STL(1) 

LDEE1 

LOA 

LOS 

AOm 

STa 

LOA 


*  0  ( 1 )  J  * 

. aores;* 

.-Ml 

Sr»l  I 
.  S  A  V  2 1 
0(  P  J 


RniiTE  JEXCHLC3  )  ilCRI 

0(1)1 

* 0 l »  s h 3 1  %      ENn  U  N  s  K  E  * 

.SAVlj*     RETUKN  TO  MAlN-PROGRAM 

S  T  C  K  J 
************************************************* 

ANGLE  FINDS  THE  TRANSFORM A T I OM-MaTrI X 
CONSISTING  OF  SINES  H  COSINES  AS  2 | 2 
DIAGONAL  MATRICES. 


.SAV1 I 

Sn3> 
24J 
*nl ' 
.specs 
0(0); 
.  an;  * 

.  A  N  T  U  N  i 
SCI  J 
NUMB  J  * 
»1  i 
SSI 
TFMPI 
NllMdl 


create  the  torn-un  pattern  for  the 
super-diagonals 


ENO  TURN-ON  PATTERN 


FIND  CORRECT  iNOEx 


0003670O 
00036800 
0003690O 
0003700O 
0003710O 
0003720" 
00037300 
0003740"! 
0003750" 
0003760" 
0003770O 
0003780" 
00037 QO" 
0003800" 
0003810  1 
0003820" 
000383oo 
0003840" 
0003850" 
0003860"> 
00038700 
0003880  0 
0003890" 
0003900" 
0003910O 
0003920" 
0003930  ) 
0003940" 
0003950" 
00039A0O 
0003970" 
0003980'^ 
0003990" 
000"000o 
0004010^ 
0004020" 
0004030O 
0004040" 
00040500 
0004060O 
0004070" 
0004080" 
0004090  0 
00041000 
00041 100 

0004120" 
0004130" 

0004140O 
0004150" 
00041600 
00041700 
00041800 
0004190" 
00042000 
0004210O 
00042200 
0004230O 
00042400 
00042500 
00042600 
0004270C 
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R'L 

cShr(  i  ) 

LDFE1 

Lt)A 

STA 

LOl(O) 
LDFE1 
LOx 
L0L(3) 

LDA 

adrn 

sTa 

RTL 
T.ShR(O) 

loeei 

LOa 
STA 

ldlc?) 

LOEEI 

LDx 

LDL(O) 

LOS 
LDL(3) 

loeei 

LOA 

rtl 

LDS 

S^rN 

sTa 
R'L 

cSmro) 

LOEEI 

LOa 

STa 


SA.lt 

11 

sc  i» 

TEMPI* 

.  ANTiJiJ  t 

TFMPI 
.  A  ')H  A  J 
*0( 3)  ; « 
*A»* 
TE^Pl S 
S  A  .  1  I 
1  > 

&roi 
«p; 
te^PI i% 

.sPEct 
*c?; 

N 1 1  m  ri  : 
. ad^a; 

.  A  N  T  U  N  ; 
Sr3; 
sst  * 

$s.-i » 

*Rt 
*<;>* 

TTMP2t 

S  A  .  1  J 
1  : 
%r  3J 

TFMP2 J* 


X 
X 

XCA 


X  A 


TEMP  CONTAINS  THL  INDEX  OF  THE  FORM 
2L*1»?L  FOR  L=0»l»2 N-l/2 


LOAn  A[?l,2l-n  FUR  ALL  I  LSS  THAN  OREQLN 
2!  U2I.2I-1  ]. 


all  pE  s  Contain  the  value  2  J  a r 2 i  • 21-11 
IN  GHC1HPS  Or  2. 


LUAD  A[  I. n  FOR  All  I 


SA  CONTAINS  AC2T-1.2I-1I  IN  EVERY  OTHER 


FOrM  AC2I-1.2I-1J-AC2I.2I] 


ALL  PE  S  OF  TEMP2  CONTAIN  A C 2 T -  1 , 2  I  -  1 ] - 

AC  2  I  ►  2  I  ]  IN  GROUPS  OF  2. 

NO"  «E  HAVE  TO  DIFFERENTIATE  BETWEEN  CAS 


**************************** 


X 

XCA 


SE  It  TEmP2=0 

LOEEI 

LOa 

LlT(O) 

IME 
NO  MASK  HUT  THE 

sEte 
setcco) 

STl(O) 

ZERE(O) 

JUMP 

LOEEI 

LlT(3) 

LDA 

CALL  SORT 

LOX 

LOL(l) 

STa 

LOa 

*********** 
SElAj  FIND  OUT 
L0L(1) 


ir2; 

TEriP2t 

=  0.0t 

SCO?* 

PART  BEYON 

I  .  A  N  0  .  E  S 

e; 

.  A  N  T  J  U  1  t  * 
AMI  it 

%rOi 
»0.5J 
tr  3t 

64(  )  ; 

Numb; 

.  AOKCt 

*o(  l)i% 

TrMPl i 
*********** 
WHTRL  tempi 

.  7.  E  H  0  i 


SET  I  WHERE  MANTISSA  IS  ZEPO» 
0  ThE  BOUNDARY  n 


IE   none  OE  THE  TLMP2  rt^RE  ZERn.  THEN  Go 
TO  ANlt  OTHERWISE  CONTINUE. 


COSINES  »   SQRT(0.5)  WHERE  TEmp2*0. 

****** 
LSS  0  .HERE  TEMP2.0. 


0004260"! 

0004290o! 

0004300^1 

0004310^! 

0004320^: 

00013300 

00043400 . 

00043500 

0004360^ 

0004370  3 

/200043800 

0004390'i 

0004400r 

0004410r 

0004420" 

0004430  i 

0004440^' 

0004450"1 

0004460^ 

0004470" 

0004480 r 

0004490r 

0004500'' 

0004510^ 

000^520^ 

0004530" 

PE00045'iO^ 

0004550": 

0004560" 

0004570' 

00045*0  - 

0004590  ■ 

0004600  ~ 

0004610  i 

0004620" 

0004630" 

00046a0 

Es0004650^ 

0004660" 

0004670" 

0004680^ 

00046901 

0004700" 

0004710  ' 

0004720" 

0004730O 

0  0  0  4  7  a  0 

0004750  " 

00047600 

0004770" 

0004780" 

0004790" 

0004800' 

0004S10" 

0004  8  2  0" 

0004830r 

0004840  1 

0004850"! 

00048*0" 

0004870" 

0004880r 
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1*1 

SETCd  ) 
CAnO( 1 ) 
ZErTU) 

LDl(2) 

CShR(2) 

CAND(2) 

LDfn 

LUS 

LDL(3) 

LOA 

chsa; 
sta 

CShl(2) 

LOs 

LOA 
STA 
CUmPCC  1  ) 

cand(O) 


*C1» 

1 1 

%rOi 
»  AN"! »  * 

.  ANT UN  J  % 
1  I 

*C2» 

.  AD*C  J 
*0( 3)  } 

*0( 3)  J 
1  i 

*C2J 
T  f  M  H  1 
*0(3  )  } 
*Of  3) » 

*r  1  !  * 


IF  TLMP1  is  N0""mERE  LSS  0  WHERE  TEmP2»0 

THEN  GO  TO  AN2. OTHERWISE 

PLACE  THE  SINES  wITH  THEIR  CORRECT  SIGNS 


TEMP  CONTAINS  INDEX  FOR  SUPER -0 1  A GON ALS 


SCASE1B: 
AN2I 


*********************** 

TLmP2=0  ANU  TEMPI  GTR  UR  EQL  0, 

z  e  p  t  ( o )  ,  a  m  l ; 

LDEEl 


LDl(3) 
LDa 
LDs 

LOi (1) 
CAnO(1  ) 

L0FE1 

CHsAI 

STA 

CShR(1 ) 

LDEE1 

«?TA 

LOl(O) 


.  ADRCJ 
*o( 3)  i 
TFMPJ 
.  A  N  T  U  in  I 

sro; 

%r\l 

<o( 3); 
i  I 

Sri  ; 
*n( 3) ; 

.  a  n  T 1 1  n  i ; 


*********************** 


XCASE2'  FiNn  CusINE 
AN1IJ     LDl(1) 

CEyOR(O) 

L0EE1 

LOA 

DVRN 

STa 

mLrN 

LlT(3) 

ADRN 

CALL 

LDS 

1. 1  T  C  3  ) 

LOA 

nVRN 

STa 

LlT(3) 

LDa 

aOrn 

c aLi_  s»RT6a( ); 

LDx        nmmbj 

ldl(3)    .aowc; 

STa        *0(3)jX 


l     SINE 
. SPEC  J 
Sri  J 
SrOJ 
TFHPl; 
TEMP2: 
TTMP1  i 
SA  : 

■1.1 

*r  3  J 
SdRT64(  )', 

sa; 

»n.bJ 
SO  3> 

*s» 

TFMP2J 
Sr3J 

ttmp?i 


Fnw  TE^2.0 


COSINES  STORED  FOR  TEMP2.0 


0004890- 

00049000 

0004910<"; 

0004920" 

0004530" 

0004940" 

0004950 

0004960' 

0004970r 

00049^0" 

0004990' 

0005000" 

0005010" 

0005020  " 

0005030" 

0005040" 

0005050" 

0005060" 

0005070 

0005050 

0005090  ' 

0005100- 

0005110 

0005120 " 

0005130  > 

0005140  i 

0005150 

0005160" 

0005170 

0005180  * 

0005190 

0005200" 

0005210 

0005220- 

0005230- 

0005240" 

0005250" 

0005260" 

0005270^ 

0005280 

0005290  ' 

0005300" 

0005310" 

0005320" 

0005330 

0005340" 

0005350-1 

0005360" 

0005370  " 

0005380  • 

0005390" 

0005400" 

0005410" 

0005420" 

0005430- 

0005440" 

000545O- 

0005460" 

0005470" 

00054800 

0005490T 
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t 

ICAS 


% 

XCAS 
AN3l 


AN4I 


LI 

LO 

SB 

CA 

**** 

E2A|  TE 
ST 
LD 
IS 
SE 
CA 
Zl 
LD 
LO 
LO 
LD 
CA 
LO 
LD 
CH 
ST 
CS 
LD 
ST 
CO 
CA 
**** 

E2tf|  TE 
LO 
LD 
LO 
LD 
CA 
LO 
LD 
ST 
CS 
LD 
CH 
ST 
LO 

E^ 

******* 


T(  3) 

A 
RN 

IL  SO 

*  *  *  *  * 
MP1   L 

A 

A 

N  I 

TC(  1  ) 

NO(  1  ) 

RTCl) 

EE1 

A 

S 

L(2) 

ND(2> 

rri 

I.C3) 

sa; 

A 

HR(2) 

EE1 

A 

MDC(  1 
NO(O) 

*  *  *  *  * 

mPI  G 

EE1 

S 

A 

L<n 

NO(l  ) 

EE1 

L(3) 

A 

hR(1) 

EE1 

SA  J 

A 

LC3) 

CHL(3 

*  *  *  *  * 


■  1.SI 

*r  3  J 

TFMP25 
RT6«(  )  5 
*********** 
SS  0 

TTMP2J 
TEMPI! 

1 1 

SrOJ 
•  a ^3  J  * 

Rr.l  S 
TFMP2; 
T  r  m  P ; 

.  ANTtjrj; 
*M  ; 

.  AO*C * 

10( 3)  J 

1  ? 
SC?» 

*o( 3 ) ; 


X  ** 

MUlTPLl   FlLLJ 

I 

« 


)  J 

*Cl JZtRT 

*********** 

Ed  o 

SrO  ; 

temp: 
TfmP2: 
,  A  M  T  U  -N  J 

?.  -  1  J 

.  4nKC> 
*0(3)  J 

l  : 

tr  i ; 

«o( 3) ; 
.«;aVi; 
)  MC«J 

*********** 


MULl  l 


STL(3) 
LlTUJ 

CAon( 1 ) 

cRnTL( 1  ) 

LDl(3) 

LOfEl 

LOy 

LOl<3> 

CAnD(3) 

LDA 

LlT(O) 

CADO(O) 

CRpTL(O) 

LOS 


.SAVl ; % 
=  r\  ,  1  ,  ;)  I 

2H  l 

.SPEC! 
*r3»* 
NllMHJ  * 

. AnKESl* 
<r\  i 
0(3  )>* 

an.  1  .0! 

?a» 
=  n; 


0005500 
0005510^ 
0005520" 
0005530 

***  000554O- 

0005550" 
0005560  ' 
0005570  ) 
0005580- 
0005590' 
0005600 
NONE  OF    ThE  TE^Pl  WERE  LSS  0»  $n  GO  TO  AN30005610T 

0005*20 
0005630" 
0005640" 
0005650^ 
0005660" 
0005670' 
0005680" 
00056O0 
0005700  ■' 
000571O  i 
0005720' 
0005730- 
0005740'. 

(0)  #AN«»  0005750 

***  0005760 

0005770 
0005780  i 
0005790" 
0005800 
0005«1 0' 
0005820  " 
0005830" 
0005840' 
0005*50  i 
0005860- 
0005870" 
0005880" 
0005890  ' 
0005900  ' 
00059m 

********************************************«no05920 

0005930- 

0005940r 

00059500 

0005960' 

0005970- 

0005980 

0005990 

0006000 

0006010 

0006020 

0006030 

0006040' 

0006050' 

0006060' 

0006070 

0006080 

0006090 

0006100' 


MULTIPLICATION  nF  T> 
SET-UP  IS  SUCH  THAT 
THE  riAlRIX  IS  TRFATEO 
SAvE  RETURN  AQORLSS 


PE-NUM4ERS 
"AMATRIX-RAsE 

LOAi)  ROW  or  AMATHIX 


0  MATslCEs.  THF 

THF  aqoress  Of 

AS  A   VARIABLE 
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MUL1H 


MUI_3i 


MuL2l 


ID|_(3) 

mLrn 

AORN 
STA 

sete 

SEtEI 

CLraj 

LOS 

LDl(3) 

LDEE1 

LOa 

LDLC3) 

£QLXE( 3) 

CAOD(O) 

CSnB(O) 

L0l(2) 

CSi|R(2) 

JUMP 

LD|_(2) 

STl(2) 

CLcC?)l 

SLITO)    i 

STA 

sete 

SETE1 
CLrA: 
LOs 

LDi_C  3) 
LDfEI 
LOa 

CLC(3)> 
SLIT(3)    i 

LDS 

LDx 

SETE 

SETE1 

CLRA  I 

LDS 

LDl(3) 

loeei 

LOs 

LOA 

TXlFM(O) 

JUMP 

L0L(2) 

CAnn(2) 

STS 

TX|_FM(  1  ) 

JUMP 

LlT(O) 

cAon(O) 

LDL(3) 
CAi)0(3) 
L0L(2) 
CAnO(2) 

LOa 
sTa 

TXfEM(O) 
LDl(3) 


BMATHlx-BASp 
MULTIPLY    BMATRIX 
FORM    PARTIAL    SUM 


.AOHESl U 
*0(3  )  1% 
*S»  % 

ma"jTj 
E,OR.-E> 
E ,  A  N  n .  E  J 


Sa  : 
.<*pEC» 

SC3J 

sb: 

Snl»MiJL3;jt       IF    *C0=1»     CHECK    IF    .SAV9«1 

.S4V9; 
*nl » ■  1 U  L  3 : 
S03! 
SOI  ? 
.  MmOJ 
Sol  5 
miil<*  : 
.  n  m  E ; 
.SAV2; 

RDi|TE;EXCHl(3)     SICrJsKIP    ,oj 
TE'-'PS 
E.OR.-EJ 

E. ANO.E5 

SA  I 
Sr3; 


ROijTE;ExCHL(3) 
Sa  t 

ss; 

E.n-».-E; 
E. AND. El 

SA  I 

. spec; 

Sr3; 

'amT  : 
T  F  M  P : 

.♦1  i 

MiiLl  i 

. aoko: 
sci  t 

0(?)J 

.♦1» 

MlJL  * 

sO» 1 »  0  J 

?4! 

.  Ar,KES2> 
tr  0  i 
. AOKOJ 

sco; 
0(?; » 

0(  3) » 
»muL2I 
.SAVl ; 


JlCRJ 


00061101 

0006120ii 

00061300. 

0006140": 

0006150"! 

0006160" 

0006170" 

0006180^ 

0006190^ 

0006200" 

0006210" 

0006220". 

0006230" 

0006240" 

0006250" 

0006260 " 

0006270" 

0006280" 

0006290" 

0006300' 

0006310" 

00063201 

0006330" 

0006340  - 

0006350" 

0006360" 

0006370" 

00063^0" 

0006390", 

0006400'"' 

0006410  ^ 

0006420  i 

0006430" 

0006440' 

0006450" 

0006460" 

0006470" 

0006480" 

0006490" 

0006500" 

0006510" 

0006520" 

0006530" 

0006540T 

00065500 

0006560,' 

0006570.-. 

0006580, 

0006590" 

0006600" 

0006610" 

0006620  ' 

0006630" 

00066401 

00066500 

0006660" 

0006670" 

0006680C 

0006690C 

0006700C 

0006710U 
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FXCHL(3)   MrHJ 


X    * 

SAMUL 
X 

*  OF 

*  pUT 
X  OF 

X  PUT 
X  AND 
X  THI 

X  MAS 
X  FXP 

%tM$t 


SAL'l l 
X   DI 
X   LO 


X  THE 
X  MUL 

SAl« 
X  TEM 


X  ADJ 

X  OF 

SALll 


X  FORM 

X 

x  con 

X  TE 


T  Hi  S  PRflC 
THE  TRANSF 
El)  IN  TuF 
THE  SYmmET 
LS  UnLy  Th 
FILLS  THE 
S  mEan$  ru 
»  T  h  U  s  •  B 
LAnATIo^S 
IfKMfllf 
SETE 
StTFl 
CLRA8 
STA 

STL(3) 
LlT(O) 
LOLd) 
CSHR( 1 

CAnoco 

CRnTL( 
LnL(2 

LUFE1 

LDX 
AGONAL  EEM 

LUL(1 ) 

LDa 
aoeo 

eqlxtc 

STi_(0) 
CLC(3) 

LOCAT JON 
TlPLlEn  RY 
SLITC3 
STA 
PUrArILy  I 
LUA 
EULXT(O) 

stl(o) 

CLC(  3)  J 
SLIT(3) 

UST  THE  IN 
THE  MATRIX 

LOS 

LDL«1) 
LOA 

MLRN 

AORN 

THE  PART 

STA 
TAKE  THE 
AL  An  MULT 
UlA'iflNAl. 

LUA 

LULC1  ) 

STL<D 

CLr( 3>» 

SLlTf 3) 


EOURF  MULTIPLIES  ANMAT(TR.1»  ThE  TRANSPOSE 
ORMA  T  I  fiN-MA  TRI  X  BY  «ASF*RASE  X  ANMAT  COM- 
pPilCFnijRE  tuLTpL.  SAMnL  TAKES  A  uv  A  gT  A^.E 
RY  Or  THE  RESULTING  MATRIX  AnH.  THUS.  COm- 
E  ELE'<ENTS  OF  THE  MRS  N/?  +  l  DIAGONALS 

REST  .-jITh  The  JUST  CALCULATEn  tLE'-'ENTS. 
LL  EFFICIENCY  .  SNCE  THE  NijMuFR  OF  MuLTIPL. 
EEN  RrnUCEO  FROM  N  X  m  TO  3(N/?*1).  FURTHER 
FOLLHw  As  THE  ALG(jRiTnM  GOES  AL"Nr, . 
ttstiia*»iH8*«t»**tt*»itttat*$»»»»t*ttt»ttMn 

F .  n  k  .  -  e  ; 

F . AND ,E  » 

TFmPJ 

.SAVlj*   SAVE  THE  RETURN  AnURESS 

=  n,l.  0',%       SET  UP  THE  LOUP  CUUNT 

.  m  ; 
)    1  :  r  sci  »  N/2 

)    hr  1  ! 
0  )      ?U  l 

)    ,<;PEC|X    SC2  =  TURN  ON  PATTERN 

Nl|r4d;<     SX  =  iNnEX  T()  PICK  (JP  THE 
LNTS  OF  A  N  m  A  T  (  T  R  ,  ) 

.A0KE51S*   AORESl=AnnRFSS  OF  ANMAT(TR.) 
♦0(1)1*    THE  DIAGUNAL  ELEMtNTS  ARE 

0  )    i  n  0  .  S  A  L  » 

.  <u  V2  : 
U  MATCH  THESE  ELEMENTS  WITH 

OF  THF  ELEMENTS  OF  BASE.  ThEy  NEED  TO  HE 

t 

)     sROilTE  SEXCHL(  3)     X  1  C  p  ; 

tfmpiix   save  the  elements  of  anmahtr.) 

n  tempi  : 
*x  : 

suo. sal l ; 

.  S  A  V  2  J 

sR0UTF.;EXCHL(3)  SICK; 

DE*  For  ThE  PROPER  PICK  UP  OF  THE  ELEMENTS 
BASE 
&A  : 
•AnRES;<   . ADRES=AonRESS  OF  mASE 
«0( 1  )J 
TEMPI  J 
TE^p  : 

1  AL  PPUOUCT. 

TEmp; 
SljPERPlACujNAL  OF  ANMAT  ARCh/E  ThE  OIA- 
IPLY  IT  MY  THE  SUPEKDI  AGONAL  OF  RASE  RELO-m 

NIIMri  ; 
.flNE» 

.  SAV-2j 

=  ROi)TFiEXCHL(  3)  *ICR» 

tAl  XAOJUST  THE  INOFX 


n006f 200| 
-0006730^! 
000674OOL 
00067501, 
0006760V 
0006770V. 
000678no' 
00O679OV 
0006800": 
000681  Oi, 
0006820^ 
0006830  ) 
0006B«Oi 
0  0  0  6  8  5  0  'i  j 
0006860'*' 
00068701 
0006880'c 
00068901: 
00069O0  1 
0006910'1 
0006920-v 
0006930  ' 
00069uo,; 
0006950  ' 
00069iSOV 
0006970"; 
0006  930  V 
0006990* 
0007000.il 
0007010  •-.' 
0007020'; 
0007030^- 
0007040V 
0007050": 
00070601, 
0007070ii 
00070801 
0007090"- 
0007100V 
00071  I0<i: 
00071201 
000M3O"- 

0007laov 

0007150V 

0007160V 

0007170V 

00071*0v 

0007190  '' 

0007200"' 

0007210"S 

0007220V 

0007230V 

0007240V 

0007250"' 

0007260" 

0007270r; 

0007280V 

0007290V 

0007300") 

0007310"! 

0007320r' 


-82- 


LUX  SSI 

L0L<1>  .AnRESlJ 

LLlA  *0(1)> 

LnL(l)   ,NHn» 

STL(l)     .SAV2JX    3RING  THESE  ELEMENTS  INTO 
X  THE  LOCTlnN  OF  THE  ELEMENTS  OF  hASE. 

C  l  c  (  3 )  ; 

SLI T{ 3 )  =R0UTe;EXchL(3)  Mc** 

E*JL*T(0)    i00»SAL2; 

STL(O)     .SAV2S 

CLC<3) J 

SLI  T  f  3  )  =Rr)UTrjLXCHL(3)  McK> 

SMH      .o: 
SAL2I   STA        TEMPI? 

EULXT(O)   SH0.SAL3; 
LOL( 1  )     Sro j 
CSUB<1)    $0l! 

JUMP       SAL"}*  FIND  THE  ROUTING  DISTANCE  FOR  THE  IN. 
%    OE  X  FR  DIAGONAL  P  T  r  K " U  P  UF  THE  ELEMENTS  OF  oASE. 
SAL3t  t-OL  C  1  5     .NMOJ 
SAL^l'    STL(l)      .SAV2! 
LDA       NUMBi 

CLC(3)> 
SLITC3)  =K0DTE;EXCHL(3)  SIcHJ 
LOS       SAJ 

ldl(1)   .aurfs; 
Lda      «o(i)i 

MLRN      TEMPI  ) 
AnKN       TFMP; 
STA        TEMPj 

Lda      m  i  )  m  r  . 

LDL(  1  )  .NMnj 

STL(l)  .SAV?: 

CLC( 3i  J 

SlIT^3)  =Rni)TE:F*CHL(  3)  JlCrfJ 

x  pick  up  The  elements  below  the  mai-m-di agonal  of  anmaktr.) 

%    AD  MULTIPLY  THEM  Hv  THE  ELEMENTS  AROVE  THE  DIAGONAL 
I    Of     H*SE 

LnS        SA; 

ldx      ss; 

LOL^l)  .AijRrsi; 

Lda  *o( n  i 

ldL(1)  sco; 

caodc  1  )  sd  1  i 

STL(l)     .SAV2J 

CLCOlji;         match  up  the  elements  with  THE 
*   EEMtNTs  of  base. 

SLIT(3)  =H0UTr;LXCHL(3)  f-ICHJ 
STA        TEMPI, 

LDA        NtJMKi 
CLC(3jjX  ADJUST  ThE  INDEX  FOR  RASE  PICK-UP 

SLlT<3)  =ROUTE,'EXCHL(  3)  *Ir;K> 
LnS        SAJ 

LDlCH     .anHES; 
LpA        #0(l)j 
MLRN      tempi  J 
ADRN      tfmpj 
I  NOW  WE  HAVE  TO  STOpF  THE  ELEMENTS  PROPERLY. 
STa        TfmPj 
L  0  a         M 1 1 M  b  * 
E«lxT(0)   $no»SAL5; 


000733O-' 

0007340  ■ 

0007350-" 

000736O" 

0007370- 

0007380' 

0007390' 

0007400" 

0007410' 

0007420- 

0007430' 

0007440^ 

0007450  • 

0007460' 

0007470  ' 

0007480 

0007490' 

0007500" 

0007510r 

0007520" 

000753O' 

0007540' 

0007550' 

0007560" 

0007570 

0007530  ' 

0007590 

0007600* 

0007610- 

0007620' 

0007630' 

0007640' 

0007650" 

0007660' 

0007670' 

0007680- 

0007690r 

0007700" 

0007710^ 

0007720 

0007730" 

0007740' 

0007750r 

0007760' 

0007770^ 

0007780- 

0007790- 

0007800^ 

0007810- 

0007820^ 

0007830'^ 

0007840'' 

0007850' 

0007860.' 

000787O*- 

0007880. 

0007890' 

00079000 

0007910'' 

0007920O 

0007930C 
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$T(_(0)     .SAV2; 
CLc(3)  J 
SLIT(3) 
SAL5I     LOs 

LD|_(1  ) 

LDa 

STa 

JUMP 
LDA 

LOK  3) 
C5u^( 3) 
STL( 3) 
CLc(3)> 
SLITC3) 
STA 

LDa 

CLC(3); 
SLITC3) 

TXET(0) 

LDS 

STA 

jump 
SAl^I     CShR(2) 

CAniO(2) 

CSHL(2) 

LDfEl 

LDs 

STa 
SAL8««    S^TE 

SEtEI 

CLpA  1 

STa        TF'<pi 

T  X  r  T  M  (  0  )   ,  +  1  ; 

JUMP       SAL'" 
X  REPLACE.  The  2k-1.?k  Awn  THF 
f  ORAKY  SoHA^iE  MATRI*  HV  FOAT 

L  I  T  (  0  )     a  y  ,  0  »  0  ; 

CRoTK(O)   ?4» 

cAnn(o)   %r>?; 
cRnTL(O)  2a; 
cLr(2); 

SALlO«    LD|.(1)     .F'JH! 

r  s  h  R ( l)   O(o); 

CExOKC 2)       %  CI  ? 
T*e;FM(0)        »SftLlOI 

ldfei  sr?; 

LUa  NJUMri; 

LlT(O)  «l ; 

A  0  m  S  C  0 ' 

lds  *a; 

CHmPC(2)J 

CAnO(2)  Sn5> 

ldeei  srtZJ 

LDa  NllMHi 

S  B  M  $  r.  o  » 

LDs  <a: 
«    REGISTER    *s    CONTAINS    THE 

LD|.(2)  .SPEC; 

LDFEt  lr2) 

LlT(O)  =0.u{ 


"KOnTL JEXCHLC 3)     MCR» 

sa; 

.aOnEs?**  aukes^*aDrfss  nr  trasE  a  temp,  matrix 

ttmp: 

*0(  1 ) ; 

%  0 0 »  +  11 

salo  ; 
n  1 1  m  ti ; 
.  n; 

*r.o; 
.saV?; 


=  R0iitL;Excml(3) 
t  r  m  h  1  ; 
t  f.  h  p  J 

srhiitl;  fxChl(  3) 
*co»  s^l  f ; 
tfmpi ; 

SAL*'  ' 

or  ci) ; 

*ns; 
0(0); 

TF^^l  '• 
*n( 1 ) ; 

F.nK.-E! 
F.  a^h.f; 


i  I C  h  i 


i  I  C  K  » 


2K  .2K-1 
ZErtiitS  . 


positions  or  temp- 


p  A  TTLPii  1  .0»  3»?»5»4» 


0007940"' 

00079501 

0007960"; 

0007970" 

0007980^ 

0007990" 

oooeooo'1 
oooftoi o  i 

00080200 
00080301 

oooeoao  -1 

OOOHOSO'i 

00080600 

000«070  ' 

00080*0  > 

0008090' 

0006100^ 

0008110"' 

0008120" 

00081300 

0008140" 

0008150  > 

0008160  i 

0008170'" 

0008180> 

0008190" 

0008200^ 

0008210" 

0008220" 

0008230" 

00082401 

0008250" 

0000260" 

0008270"* 

0008280' 

0008290'. 

00083001 

0008310  i 

0008320" 

0008330' 

0  0  0  8  3  a  0  '■ 

0008350" 

00083601 

0008370V 

0008380  V 

00083901 

0008400^ 

0008410.- 

0008420" 

0008430" 

0008440" 

0008450" 

000  8460" 

00084701 

0008480" 

0008490O 

0008500T 

00085101 

0008520" 

0008530i' 

0008540O 
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STA 
*  THE  MATpIX  HAS 
X  MATKlX  STnREO  I 
LlT(O) 
CRnTR(O) 
CAnn(0) 

CRnTL(O) 
SAllU    L0L(1) 
L0L(3) 

cAnoO) 

loa 

CAnn( 1) 
STa 
TXpr'UO) 

LOL(3) 

EXrHL(3) 

X####**#*###»tf *** 


A0DIT« 

X 

X 

X 


FILL! 


AOI» 


ADIH  « 


.AOKLS2J 

eo(  1  >  J 
mow  its  final  Form.  *e  now  map  THE 
N  TsASE  INTO  RAsE. 

=  1  .  0.  ii ; 

?4j 

2  a ; 

.  aoP.ESJ 

.  4nRE'S2  J 

Of  3)  J 
SCO? 
0(1  )J 
•sali l ; 

.SAV1  ; 

STr«J 

»tneuaa.tvntmtttti%tt»tt«tt*nttttttttttt 

ADnIT  rALCOLATES  THE  SUM  Or  The  OFF- 
OlAGUNALS  SiUARE  ANO  DIVIDES  TT  BY  THE 
SUM  UF  THE  DIAGONALS  SQUARE. 


STL(3) 
LlT(l) 
STl(1  ) 
LlT(O) 

cAonco) 

CRflTL(O) 

LDLC3) 

StTE 

SETFl 

CLrAI 

LOS 

ldefi 

LPl.(2) 

CApnc?) 

IDA 

mLrN 
AOrN 

LDS 

TXrFM(O) 

SETE 

SEtEI 

LlT(O) 

LDS 
RTL 

LDS 

AORM 

CAnn(o) 

LlT(l) 
EOLXF(O) 

ldlci) 

FQLXTC 1 ) 

STa 
CLpaj 

LDs 

LDLd) 

LDEE1 

LDX 

LDl(2) 


.  s  a  v  i  ; 
=  o»s 
.SAV5; 
=  o.  1 .0;  t, 
Sn3! 

?a  > 
.spec; 

E.nR.-E? 
E.  ANM.f. # 


CHEcK  FOR  nFF-OTAGONAL  OR  nlAOONAL  SUM 
ADDING  THE  ROWS  INCLUDING  A  C  I  •  I  3 • 


SA  I 

*r3 

.  An 
SCO 
0(2 
*A  I 
«S» 

r-A  i 

•  A  n 
F.n 
E.A 
=  1  » 
IM 
IS. 
«.p  : 

*S» 
«ro 
=  64 

sri 

.  SA 
Sol 
TFM 


ra; 
; 
) ; 


I i  *     ENn  ro^-sum 

R.-EJ 

ND.EJ 


*        SUM  UP  THE  ELFmFNTS  IN  EACH  PF. 
0(0);*   USED  FOR  ROTH  OFF-OUGONALS&  OIAgOnALS 


»AoI1;«  ENo  LOr,-SUM 
V5J 

» Aui2; 
p; 


SAS 

.spec; < 
f  c  l » 

nomh  ; 
. adha; 


PICK  UP  THE  A [ I • I ]  S, 


00085500 
00085600 

00085700 
0008580" 
0008590" 
00086000 
0008610" 
0008620O 
0008630" 
00086400 
0008650" 
0008660" 
0008670" 
0008680" 
0O0B69OD 
0008700" 
00087100 
0008720" 
0008730" 
0008740" 
0008750" 
00087600 
000377OO 

0008780" 
00087900 

000 8800  0 
0008810" 
0008820" 
0008830" 
0008840" 
0008850" 
0008860" 
0008870 o 
00088*0" 
0008890" 
0008900" 
0008910" 
0008920^ 
0008930" 
0008940'1 
0008950) 
0008960^ 
0008970" 
00089H0" 
0008990^ 
0009000" 
000  901 00 
0009020" 
0009030" 
0009040" 
0009050" 
0009060" 
0009070" 
0009080C 
0009090" 
00091000 
0009110" 
0009120" 
00091300 
00091400 
00091500 
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LOA 
MLR 
SET 
SET 
LIT 

L0L 
ALi 

STL 
JUM 
A  D  I  2  »  L  0  s 
LDa 
SBr 
LIT 
IMF 
SET 
ONE 
qVr 
STa 
SLl 

loa 
stl 

A0I3I     LDL 

F*C 
*   M******** 

TRASPUSi  fIl 

* 

X 

X 


u 

E 

El 

(0) 

(  1  ) 

T<1> 

(1  ) 

P 


N 
(0) 

r(0) 

ST(O) 
N 

T(0) 

0(0) 

(0) 

(  3) 

HL(3) 
*  *  *  *  * 

Li 


*n 
*A 
E, 
E. 
■  1 
.  s 
=  1 

•  S 

An 

*A 
TF 
*S 

=  n 

sr 
II 

•  A 

*s 

TF 
=  T 
if. 
.  <! 
.  9 

SI 

***  * 


(2)) 
I 

ok.-e» 

ANO.E J 

1 

A  V5  I 

II I 

I 

MP) 

J* 
,o» 

OJ 

n!3< 
:  * 

V1PI 

EMP5 
OJ 

i 

AVI  I 

C*» 


SUBTHACT  THE  AtTMI. 


CONVERGENCE  F4CTUR  KSI  FOuNO« 


STL 
LDL 
LOF 
LDA 
LOS 
ADm 
STA 
CSh 

rtl 

LOE 
LDA 
STA 
LDL 
LOE 
LOS 
LOL 
LDa 

CHS 
STA 

LOL 

********* 


(3) 

(0) 

El 


RCO) 
El 

(0) 
El 

(3) 

At 

(  3)  . 

*  *  *  *  * 


r************************************** 

TrtASpOS  FlNnS  THE  TpAN$POSF  Or  THE  ANGLE' 
MATHIx  AY  JUST  CHANGING  THE  SIGNS  OF  THE 
SUPEK-MAGHNAL  ELEMENTS. 

. SA Vl  J 

.4nTij.m;J    FInj  THE  IndEX  FUR  ABOVE  Elements 

SCO  I 

Nl|Mrf  I 

=  1  I 
«S» 

TE^'LOA  NUM«» 
1  I 
84.  1  I 

UO! 

tfmpj 

.SPECi 

vo; 
Tempi 
•  a  n  H  C  i 
«n(3) i 

«n(i) t 

SAVl JEXChL( 3)  MC*> 


*  I.NnEX     IS    FOUNT 

i 


**** 


CONVi 

1! 

X 


FlLL» 


************************************************** 


CONl 


STi C  3) 

LlT(O) 

cAnoco ) 

CSu°(0) 
CRnTL(O) 

CLcd)> 
L0L(2) 
CShR(2) 
cExPRC 1 ) 


.SAVU 
»n. 2. o  J 
So-}  J 
*ni  ; 
24  J 

,  F  «i  H I 
0(0)  I 

8C?» 


CONV  CHECKS  FOR  THE  A8S. -VALUE  OF  THE  OFF 
DlAbUNALS  ELEMEMS  LSS  .ftD(INopx). 


0009160."! 

0009170"! 

0009180"! 

0009190' ' 

0009200"! 

00092100. 

0009220": 

0009230"* 

0009240" 

0009250" 

0009260". 

0009270  i 

0009280" 

0009290" 

0009300" 

0009310i' 

0009320' 

0009330" 

0009340  ' 

00093501 

0009360" 

0009370  1 

0009380" 

0009390" 

0009400" 

00094100 

0009420" 

0009430" 

0009440" 

0009450"; 

0009460) 

0009470 

0009480" 

0009490" 

0009500" 

0009510 

0009520- 

0009530  i 

0009540 

0009550" 

00095*0" 

0009570" 

0009580" 

0009590" 

0009600" 

0009610" 

0009620- 

0009630" 

'0009640" 
0009650  ' 

-0009660':. 
00  09  6  7  0", 
0009680" 
0009690" 
0009700-' 
0009710- 
0009720" 
0009730- 
0009740" 
0009750" 
0009760. 
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TXfFM(O) 
RTL 

csmr( 1) 

IDEE1 

LDs 

CSwL(l) 

LOEF1 

L0L(2) 

ADM 

LOS 

LOLH) 

LDfEl 

LDL(2) 

LOa 

LDl(O) 

Ml 

StTE 

SETE1 

CHSA» 

COMPCd  ) 

LDpEl 

CLRAJ 

CUmPCC 1 ) 

LOfEl 

LOl(O) 

L0l(2) 

ial 

SETCC2) 
STl(2) 
LDlO) 
EXcHL( 3) 
************** 


•  CONJ 
Scl  I 

NIJMBJ 

$A»n 

it 
sen 

SPJ 

1 J 

SCU 

.ONE  J 
SC2> 

*  A  I 

.SPEC  J 

set  I 

.  ADKAJ 
*0(2) J 

.zeho; 
sco; 

I  .  A  N  1 .  Z  I 

E ,  a  n  q  .  e  ; 


Sri  ; 


fill; 


TRPSI 

% 

t      ORIGINAL 

stl 

LIT 

cad 
cRn 
ldl 
loe 
lDa 
e°l 

STL 
CLC 
SLI 
LDS 
LDX 
LDL 

LDa 
eol 

LDL 
CSii 
STL 
CLC 
SLI 
STA 
LOA 


TS3I  t 


TSi 


MATH  I 

(3) 

titit 

(0) 

0(0) 

TLCO) 

(3) 

El 

XT(O) 
(0) 
(3)  J 
T(3) 


(3) 

XTCO) 

(3) 

B(  3) 

(3) 

(3)1 

T(3) 


TSll 


«ci  > 

,SAVM 

. jn( o)  i 

V25 
II 

.rnNVE; 
.SAVi ; 
stch; 

A************************************** 

THIS  ^ROCEOURE  TRAMSPOSES  THE 
X. 

.SAVl;X   SAVE  KETURN  AORESS 

»f).l»Ol<   sET  UP  L00p-I.M0EX  I 

2a; 

.SPEC! 

Sr3l 

Ml|Mtf  I 

iSnOiTS! 

,saV2;*  riute  o  »  i 

=  R0i)TE;EXCHL(3  )  JlCRJ 
*A» 

ssi 

.flORB!?   4ORESS  OF  eigv 
*0(3)l 

$to»  tsi> 
*co* 

.5AV2I*    ROUTE  D«  N-l 

=  R0i!TE;EXCHL(3)    UCRJ 
TEMPJ 
N11MBI 


************** 


00 
00 
00 
00 
00 
00 
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APPENDIX  C 
The   Subroutine  EEERL,   Jacobi-Like  Method 
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CAPnf ?■> 
L  I  T  (  3  ) 
FQLyF(9> 
*  SS  AND  SA  rrwTAT 
L  I  T  f  0  ) 

ime 

SFTrfO) 
PNrSFf 0) 
JUMP 
LDFrl 
LDA 
* 

IAL 
SFTr 
SFTrfO) 
I  FAnPfOW 
L  TTf ^) 
C  A  N  P  (  0  5 
STLrO) 
SL  iTfl  ) 
OADnf 1 ) 
LPAP( 1 ) 
STL  <■  ^5 
CLCM  51 
SI  TTf  1  ) 
CADPf 1 > 
LPAPf 1 ) 
STl  r3) 
FND3t    LOL^ 

FXfMl  (  3") 
% tittttt t it*  *  tit* t 
ML  TRPS  J   Fill! 
I! 
t 
t 
t 

STi  f3) 
I  IT<0> 
CAnP(0-> 
CROTLfO) 
MTl  J  t      SFTP 
SFTF1 
CI RA|* 
LD* 
STA 

I  Dl  (?) 
LDFF1 
I  IT(?S 
CAPPC?} 
CRPTLf?) 
MT!        LDA 


-II 

ss>o(2)i 

SRI 
SSI 

I.AI 
F.AI 
SS» 
F.OI 
F.AI 
SAI 
S02; 
=  64; 

tC3,FMP2jx     THF  LARGEST  TS  FOUND 
N  IT. 
=  0.1 

SCO;?!    PHFCK  FOR  THF  LARGFST  FLFM*.  TO  BF  ZFRD 
II 

fnd: 
$ci 

TFMI 


i  N  P  .  F I 

iNP'.Fl 

IR.-FI 

iMO.ri 

'i 
i 


i 

31 

I 


<PI*     FlNP  ROW-  AMD  CPIUMW-imDEX  RY 
COMPARING  TFMP  AGAINST  $S. 
SSI 

-I .ANP.FI 
Fl 


=77i8j 

SC3I 

.  M  A  X  R  I  ? 

=TFMP) 

SCOJ 

SC3I 

,GI* 


PCWINPFX  FOUNP 


STPPF  THF  LARr-FST  ELEmFNT 


=TFMP1 I 
SCOI 
SC3I 

.MAXCl*        THP  CPLUMNINPr  FPR  ABOVE  FLFM. 
.  S  A  V 1  I 
SICRI 
tttt1*1t11ttt1t1tftt111**t111t**tt111t 


MULTIPLY  TwO  MATRTCFS,  ONr  OF 
wHICH  IS  THE  TRANCPOSF  OF  THF 
OTHFP,  WITHOUT  ACTUALLY  TRANS- 
POSING THF  MATRIyI 
SAVE  RFTURN-APRESS 
SFT  UP  LOPP  FOR  PPW-FrTCH  OF  PA^F 


•  SAV1 

~Op\p 

SP3I 

24J 

F.OR. 

F.  AND 

SAI 

TEMPI 
.SPFC 
SC2I 
=  0#1# 

SD3J* 
24j* 
NUMB  I 


01? 


'.Fl 


CLEAN  VP    SA  *  SS 


01 


TNMFP  LPOP  FOR  PICK-HP  OF  TRANS- 
POSE FLFMFNTS  OF  BASF 
*    MAKF  FPPPFR  TNPEXTNG 


00094600 
00024700 
00024800 
00024900 
00025000 
00025100 
00095200 
00025300 
00025400 
00025500 
00025600 
00025700 
00025800 
00095900 
00096000 
00096100 
00026200 
00026300 
00026400 
OP026500 
00026600 
00026700 
00096800 
00026900 
00027000 
00027100 
00027200 
00027300 
00097400 
00027S00 
00027600 
00027700 
00027800 
00027900 
O0028000 
00028100 
00028200 
00028300 
00028400 
00028500 
00098*00 
00028700 
00028800 
00028900 
00029000 
00029100 
00099200 
00029300 
00029400 
00029500 
00029600 
00029700 
00029800 
00029900 
00030000 
00030100 
00030200 
00030300 
00030400 
00030500 
00030600 
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%t*1 

TRPS 

* 

X 

X 

X 

X 

x#c* 
X 


TS31 


TS» 


STI  f?). 
Clrf 3)1 
SLtT(3) 

LD* 
LDY 

LDl  f  1) 
CAnPM  "> 
LOa 

LDi  (1) 
MLRM 
ADRN 
Clr(3)» 
LDl  (I) 
STI  (  1  ) 
SL  TT(3> 
STA 

TXFFM(?) 
LDA 
LDl  (3) 

txftmco) 

JUMP 
I  Dl  (3) 
FXrHL' 3) 

tlttltttft*** 
I      F  1 1  I  I 


.SAV?I 

sROUTFlFXCMLf 3) 
SA| 
SSI 

.ADRFSI 
SCOlX 
Of  1>J 
.ADRFSI 
♦0(1  )>* 
TFMPJX 


SICRI 


PlXK  UP  ROW  T  =  SC0  OF  PASF 


MULTIPLY  PY 
AnP  PARTIAL 


TPANSP. 
PRP0UCT 


f|  FM.  OF  BASF 


TSl  I 


STl  (3) 

ttitttmttft 

I    lT(0) 

r AnO(0> 

rRpUCO) 

LDl  f3) 

LDFF1 

LDa 

FOl  XTfO) 

STl  f0) 

C  L  r  C  3  )  I 

SLTTC3> 

L0< 

LDY 

LDl  (3) 

LDa 

FOl  XTCO) 

I  Dl  f  3) 

C  SliP  C  3  > 

STl  f  3  1 

r  L  r  (  3  )  I 

SLtT(3) 

sta 

LDa 

FCl  YTf05 

CI  rC3)l 

SLTTC31 

SKIP 


.ONFl  ' 
,SAV2» 
=RPUTFlFXCHl  f  3)     SICRJ 
TFMPI 
#MTJ 
TFMPI 
.ADRFSU 

sroi 

0(3)1 
#♦11 
MTll 
.SAVl » 
SlORl 
*tt1ttttt1tt*tt1111ttttttttt1ittt*11tt1t1 

THIS  PRPCFnilPF  TRAMSPOSrS 
THF  PRTfiTMAL  MATRTv*  TP  PRFPARF 
TT  FPR  TMF  MUI  TIP!.  TCATION 
(A-)TXA=AfTRAMSP.)Y(A(TRP.)TPP,'i 
AfTRP.)  IS  TE>PLY.STPPFD  IN  AumAT 

.SaVIJ*  savf  rfturm  aorfss 
ittttifttt 


=0* 1 #01 

SD3I 

?4| 

.SPFCJ 

SC3> 

NUMRI 

SPO.TSI 

.SAV2JX 


sft  vp   lppp-inpfy  t 


RPUTF  P  =  I| 


=R0UTFIFXCHL(3)  SICRI 
SAI 
SSI 

.ADRAjf  APRFSS  OF  PASF 
*0(3)l 
SDCTSU 
.Ml 
SCO! 
.SAV2I*  RPUTF  D»^-T 

=  P0llTFIFXrML(3)  SICRI 
TFMPI 
NUMRJ 
SDO#TS?) 

=  R0UTFIFXCMI  (3)  SICRI 
#01 


00030700 
00030800 
00030900 
00031 000 
H0031100 
00031200 
00031300 
00031 400 
00031500 
00031600 
00031 700 
00031 800 
00031 900 
00032000 
00032100 
00032200 
00032300 
00032400 
O0O32S00 
O0032600 
00037700 
00032600 
00032900 
00033000 
00033100 
00033200 
O0033300 
00033400 
00033500 
00033600 
OOO337OO 
00033800 
00033900 
00034000 
00034100 
O0034200 
00034300 
000  34  4  00 
00034500 
00034600 
00034700 
O0034800 
00034900 
O0035000 
00035100 
00035200 
O0035300 
00035400 
00035500 
00035600 
O0035700 
00035800 
00035900 
00036000 
O0036100 
00036200 
00036300 
00036400 
00036500 
00036600 
00036700 
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TS?l       LD$ 

LDA 

LDI  C35 
STA 

TXFFMf 0) 
LPl  C35 
FXrHLC35 

%i*tiitiiiiiiitit 

SHUFi     FIl  Li 


STl 
%11111111111 
*  PAPT  0NF» 
LDI 
LPr 
LDI 


FOI 

JUM 

LPl 

CSli 

STl 

I  IT 

CAD 

CRP 

LDI 


PARI  : 


PARA1  I 


PARA2I 


r  AH 
LPfi 

CLr 

SI  T 
STA 
TXr 
M 
F|  Fm 
I  IT 
f  An 
CSn 
rRn 
LDi 
CAn 
LDa 
LDi 
CAn 
STA 
TXF 
I  IT 

CAn 

CPn 

CAn 

I  Di 

CAn 

LDa 

CSit 

STA 


f  3> 
11*  11 

(  3) 

Fl 

C05 

XFC05 
P 

f  15 
P(l  1 
C15 
C05 
H(Oi 
TLf  P5 
(M 

P(1  > 

f  3  5 1  y 

T(3)  i 

FMC05 
HW  RF 
FMT  I 
(0) 
P(O-) 

p(oi 

TLC05 

ri  > 
nci  i 

f  1 5 

PCM 

fmco5 

f05 

nccn 

TLCO) 
PCO> 
fl5 
PC1  5 

PCI  5 


SAJ 

TFMPJ 

.ADROlX   APRFSS  OF  TBASF 

#0C3)l 

>TS3I 

.SAVlir   RFTURN 

SICRI 

it  mm  minimi  tut  in  inn  m  tint 

THTS  PRPCFPURF  FALLS  INTO  TWP  PAPTS 
PART  OK'F  PRTNCS  RPW  MAyR  IMTP  THF  PI  ACF 
OF  THF  FIRST  ROW.  THF  SFCONn  PART 
PRIMGS  THF  CPLHMW  MAXC-MAXP  TNTP  THF 
PI  ACF  nr  THF  SFCPND  COLUMN'.  THIS 
PRPCFDMRF  IS  FNTFRFn  PNI  v  IF  MAXR+MAXC 

NFO|  r. 

.SAVll*    SAVF  RFTHRM  APRFSS 

t  timtt  tut  ti  111111111111  tint 


.SPFC? 
SC3I 

.MAXR)* 

sdo,+i;x 

PAR2J 
.NJ* 
SCOJ 
.SAV?!f 
=0* 1 >0J 
SD3* 
24; 
.APRFSJ * 

SC01X 

ocn; 


TF  HAXR=0#  THFN  THE  I  ARfiFST 
FLF^FK'T  IS  Al  PFAPv  TM  RPW  0. 
WF  JUST  NFFP  TO  AP.IMST  THF 

COLUMNS 

POHTF  PTSTANCF  n=M-MAxR 


CONTAINS  THF  ARrSS  PF  THF 
MATRIX  UNPFR  C ONS I PFRAT I PN 
PPUTF  ALL  PPWS  PF  MATRIX 


HP  TP  THF  ROUTF  PPPCFPHRF 

=routf;fxchl (3)  SICRI 

PCI  5  J 

*PAPl J 
ARRANfiF  THr  MATRIX  SP  THAT  THF  lARf.FST 
S  IN  THF  FTRST  RPW. 
=0*1*0) 

$015;*    I  I M  IT  TS  maxp-1*  sTORFD  IK'  $nl5 
*DlJ 
24; 
.  APRFSl 

sco; 

0C15) 

.ADRP»* 

SCO; 

0C1  WJ 

>PARA1 ) 

=0* 1*0JX 

SD3; 

24; 

SD15J 

.ADRFSI 

SCO; 

ocn; 

SP15) 
0C15I 


APRFSS  PF  TPASF 

STOPF  RPWS  TFMPPPARILY  IN  TRASF 
LOOP  TP  FFTCH  RPuS  MAXR  TO  N-l 


00036800 
00036900 
00037000 
00037100 
00037200 
00037300 
00037400 
00037500 
00037600 
00037700 
00037800 
00037900 
00038000 
00038100 
00038200 
00038300 
O0O38400 
00038500 
00038600 
00038700 
00038800 
00038900 
O0039000 
00039100 
00039200 
00039300 
00039400 
00039500 
00039600 
00039700 
00039800 
00039900 
00040000 
00040100 
00040200 
00040300 
00040400 
00040500 
00040600 
00040700 
00040800 
00040900 
00041000 
00041 100 
00041200 
00041 300 
00041 400 
00041500 
00041600 
00041700 
00041800 
00041900 
00042000 
0004P100 
00042200 
00042300 
00042400 
00042500 
00042600 
00042700 
00042800 
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P  A  R  A  3  t 


Itttttt 
%     PAH 
PAP?  t  t 


%      If  T 
*   pnsi 

%         IK   P 


PAPT9: 


I  PL 
CSn 

LDl. 

r  An 
LPA 
LPl 
f  An 
CAT 
STA 
TXr 
ttttt 
TWpi 
LPl 
CSn 
FOl 
Jt'V 
HF  Tr 
TlPN 
APT  P 
ST| 
I  Dl 
C  An 

I  PA 

STA 

I  IT 

f  An 

CSn 

t  Dl 

CAn 

LPA 

STa 

CSn 

TXf 

LDa 

STa 


(  3) 
P(3l 
f  1  } 

nd  ■> 
n ) 

Pd  i 

r  m  f  o ") 
tint* 

fOl 
P(O-) 

yp  fo) 

p 

ST  TS 

or  TM 
ttr\ 

(Of 

r  1  ) 

PC  1  i 


(?1 

p(?) 

P(?1 

PC  31 


p(1  > 

r  M  1 7 i 


lppp  fftchfs  thf  rows  or 

n  STPRFS  THFM  PACK  TNTP 

x  unptp  cdnsio. 


SC3  CPMTATNS  RDW  TNPFX  FOP  BASF 


TXrFKTO)   >PARA?J 
LlTCO>      *0»l»OJ 
CAnP(O)    SD15J 
CSnP(O)    $011 
CPnTLfO)  7k\%  THIS 

TPASF  AN 
THF  MATPI 

$ni 5  j  y 
# APPP t 
$co» 

oihi 

.AORFSJ 
$C3J 

sroi 
ocm 

»PARA3I 
tttttttttttttttttttttittttttttttittttitt 


. MAXC|*    FH'P 
$P15»*     WFNT 
SHI »+l J 
PVFP1 J 
SATTSFTFD.  THfN 
F  SFCPNP  CPI  UMN* 


OCT  WHFRF  THE  COLUMN 
TO  AFTTR  APPVF  PFARp". 


.SAv9?r 

. adpfs;* 
sco; 

0(1  M 

tfmpj 
=-i,i»o; 

$C0* 

$01  \ 

.  ADPFSJX 

$C?J* 

0(311 

0(1  M 

$011 

>PART?J 

TFMP: 

0(1  ")!* 


THF  COI  IJmn  TS  TM  THF_ 
RFCAUSF  OF  THF  PFARr' 

SAVF  THF  NiFw  COLUMN  INDFy 
FfTCH  POW  MAXC-M4XR 


PULL 
MAXP  P 


ROWS  APOVF  ROW  MAXC- 
PWN  RY  ONF 


pOW  MA 
rr    THF 


XC-MAXP  IS  NOW  IN  PLACF 
SFC"WO  ROW. 


%$itttt 

%     PFARR 


PART^l 


ttt* 

AKIGF. 

I  I 

rA 

CP 
SF 
SF 
CL 
LO 
LD 
I  0 
LPL 
CA 
I  0 
ST 
CL 
SL 
ST 


ttttt* 

TH*"  C 
T(05 
PP(Ol 
CTL(P') 
TF 
TF1 
PA  I 

I  (3l 
rrl 
r1  > 
nod  i 

A 

i  (0) 
r(2M 
TT  C  3-> 

A 


tttt*«*t 

Pl.llMK'S  NFXT'.PllT  SKFW  FIRST  FOR  PFTTFP  ACCFSS 
=0»1#0) 
$031 
?Ul 
F.OR.-FI 
F.ANP'.FJ 

$A» 

.SPFCl 

$C3> 

. AnPFS; 

$001 

fldll 

.SAV2J 


=routf;fxchi  (  3) 

0(1)1 


$TCRI 


00042900 
00043000 
000431 00 
00043200 
00043300 
00043400 
O0043500 
00043600 
000 0  37  00 

00043800 
00043900 
00044000 
000441 00 
O0044?00 
00044300 
00044400 
n0044500 
00044600 
00044700 
00044600 
00044900 
O00A5000 
00045100 
00045200 
00045300 
000^5400 
00045500 
OC045600 
00045700 
00045800 
00045900 
00046000 
00046100 
00046200 
00046300 
00046400 
00046500 
00046600 
00046700 
00046600 
00046900 
00047000 
00047100 
00047200 
00047300 
00047400 
00047500 
00047600 
00047700 
00047P00 
00047900 
00048000 
000481 00 
00048200 
00048300 
O0048400 
00048500 
00046600 
00048700 
00046800 
00048900 
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txffm(o) 

LDl  CO) 
STI  (0) 
I  iTfO) 

rAoO(O) 
CSHP(O) 
I  Da 

CLff 3^1 
SLtT(3) 
LDs 
I  Dv 

L ni  (l) 

IDA 

STA 
PAPT4  t    LDl  f  2) 
STl  (2) 
LDA 

CLrf 3)  J* 
SLTTC3)  = 

LDS 
STs 

LDA 

LDl  (?) 

STl  (2) 

C  L  r  (  3  )  I 

SLtT(3)  = 

STA 

LDV 

TXFFM(O) 

LDA 

LDl  f2> 

CSmPC?) 

CAr>n(?-) 

STI  f  2  > 

C  L  r  C  3  )  I 

SL  T  T  C  3  1  = 

STA 

I  It(0) 

CAnD(O) 

LDl  f  n 

FOl  XT(O) 

PARTAJ    STI (0) 
LDA 

CLr(3)t 
SLtT(3)  = 
S7a 

PART">:    ALTT(1l 
T  X  F  F  M  f  0  ) 

0VFR1  I  J   LDl  (3) 

FXrHLf 1) 

MHLTPL  I   Fli  I  J 

X 

X 

X 

STi'f  3) 
LItCI) 
CAnDd  ) 
CRnTLf 1  ) 
LDl  (3) 
LDFF1 


»PART3J 

.SAV9JX     MAXC-MAXR 

.SAV?J 

=-l#l>0) 

$P30j 

*D1  ) 

NUMRJT   ADJUST  INDEX  FOR  COL.-FfTCH 

sROUTFJEyCHLO)  HCRJ 
$AJ 

SSJ 

. AORFS* 

*Ctnj*    FETCH  COLUMN  f(MAXO-MAXP) 

TEMPI 
,NMOI 
.SAV?J 

*X»X     ADJUST  PICK-UP  IMTEX 
FDP  RPWS  .SAV9-1  TO  ! 
R0UTFJFXCHL(3)  SIfRI 
SA  X 

TFMP1 I 
f Of  1  )J 

.one» 

.SAV?I 

RPUTFJFYCHI.  (3)  SICP? 
*0f  1  )  J 


TFMP1 I 

>PARTAJ 

TFMP JX 

.N) 

S030I 

*D1  I 

.SAV?J 


RPUTFJFXCHI  f3)  STCRI 
*0  f  1  )  J 

= - 1 » l » o  i 

$D?J 
•AHPFSI 

SD2/PART5J 

0(n; 


FPU  MAXC-MAXR  TN  TO  THF  SFCOND 


UNSKFW  MATRTV 


roiitf;fxchl(3)  stcri 

0(1)1 
si) 

>PAPT6) 

.SAVlJX    RETURN  TO  MAIN  PRpGRAM 
SICRI 
t*tt1tt1t1tttt*tttt**Ut*911t1t1tt* 

MUI  TTPLICATION  pF    TVD    MATPTCFs'.  THF 
SET-UP  TS  SUCH  THAT  THF  ADDPFSS  OF 


.SAVl  It 

=0*1#0I 

SD3I 

24) 

.SPFCI 

SC3) 


THF  MATRIX  IS  TRFATFO  AS  A 
SAVF  RETURN  ADDRESS 


VARIAPLE 


00049000 
00049100 
00049200 
00049300 
00049400 
00049500 
00049600 
00049700 
00049800 
00049900 
00050000 
O0050100 
00050200 
00050300 
00050400 
00050500 
000^0600 
O0050700 
00050B00 
00050900 
00051000 
00051 100 
00051200 
00051 300 
00051400 
00051500 
00051600 
00051 700 
00051600 
00051 900 
0005POOO 
00052100 
O0052200 
00052300 
00052400 
00052500 
00052600 
00052700 
00052800 
00052900 
00053000 
000S3100 
00053200 
00053300 
00053400 
00053500 
O0OS36O0 
00053700 
00053800 
00053900 
00054000 
00054100 
00054200 
00054300 
00054400 
00054500 
00054600 
00054700 
00054800 
00054900 
00055000 
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MULl  I 


MUL1  l  » 


KUL3: 
MULflt 


I  Dx 

I  01  (  3) 

CAnPf  3) 

LDa 

I   TTf  0) 

c.  Annr  fit 

CRPTLf 0) 

LD< 

I  PI   (3*> 
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ADPK' 
STA 
SFTF 
SFTF1 
CLpA  I 
LDS 

LDI  (31 
I  DrFl 
I  Da 

rci  vf  f  oi 

LOl  (3) 
r  e  i  y  f  ( 3 1 
rAnPf o> 

CSliP  CO-) 
CSnP(Oi 
L  Oi  (  ?} 
CSliP(?1 
SKTP 
I  01  (  ?1 
STI  (?1 

CLr(3i t 

SLtT(3i    ■ 

STA 

SFTF 

SFTFl 

CLRAI 

LD< 

LDI  (3> 

LDFF1 

LDA 

C I  C  (  3  >  I 
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LDV 

SFTF 
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CLRAI 

I  OS 
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JUMP 
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0(3) 
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F.OR 
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.SAV 
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FSJf 

If 
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.-rJ 
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PF-MIHBFPS 
xAMATRIX-PASE 

LOAD  POW  OF  AMATRjy 


PMATPTX-RASF 
MULTIPLY  PMATRTV 
FOPM  PARTIAL  SUM 
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SC3! 

SXI 

ROUTFIFyCHI  (3) 
SAJ 

SSI 

F.OR.7FJ 

F.ANO'.  F) 

*A  J 

.SPFCI 

SC3I 

MANTI 

TFMPI 

,♦1  I 

mhli  1 

. APPOI 

*C\l 

0(2)1 

*  +  \  t 

MULl 

=1#0#0J 

?oj 


SICRI 


SICRI 


00055100 
00055POO 
00055300 
00055400 
00055500 
00055600 
00055700 
00055800 
00055900 
000«s6000 
000561 00 
O0056P00 
00056300 
00056400 
00056500 
00056600 
00056700 
00056800 
00056900 
00057000 
00057100 
00057?00 
O00S7300 
00057400 
00057500 
O0057600 
00057700 
00057800 
00057900 
O0058000 
O0058100 
O0058200 
000S8300 
00058400 
00058500 
00058600 
00058700 
00058800 
00058900 
00059000 
O00S9100 
000^9200 
00059300 
00059000 
00059500 
00059600 
O0059700 
00059800 
00059900 
00060000 
00060100 
O0060P00 
O0060300 
00060400 
00060500 
00060600 
00060700 
00060800 
00060900 
O0061000 
00061 100 
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MUL2I 


t  **** 

MULSA I 

X 

IF  ANM 
VAI  UES 
T  1=0*2 
FORM  T 


CAPP(0> 

CRnTL(O) 

I  Di  f  3) 

CArn(3^ 

LPi  (?) 

CAnnc?^ 

LDi 

STA 

T  X  r  F  M  f  o ) 

LDi  C3> 

FXrWL  f 3) 


?X?  MA 

so  tha 

FUPTHE 


ritti*tt 


SD3I 
24| 

.AnRrs2i 
sco; 

.AnROl 

scoi 

0(2W 
Of  3")  I 

>mui  2; 

.  SAVl  J 
MCRJ 
*********+**+*•*****+*** 

F  1 1  I  I 

MULTJPIY  AMMAT(TP.)  X  (BASF  V 
AT(TR.>  CONSISTS  OF  COSH  ft  SlNH*  WE  USE  THESF  TWO 

PK'J  Y  Tp  MULTIPLY  ROW  I  AND  ROW  T  +  l*  WHFRE 
Pb»\\  .*N-?.  THFM  WE  ADD  ROW  I  ft  ROW  1  +  1  TO 
HF  MFW  ROW  T  OF  PASF. 

Ir  ANMATfTR.)  CONSISTS  OF  COS  ft  SIm,  TMFN  TMF 
TPIrFS*  GOING  DOWN  TMF  DIAGONAL*  WILL  BF  DlFFFRENT 
T  Tl-'F  CDS  *  STN  MAVF  TO  PF  Pill  LED  AS  I  INCRFASES'. 
R  FvPIAmATIOkS  Jm  TMF  PROGRAM. 
1tt1tt*tttt11l1ttt1it*1ttttf*1ttt11tt1itt*1*9ttfH$ 

r  T  ■    r  1  s  c   a  \i  1    i  a  CAl/r      Tnr      nfTnDii       «  n  n  r>  r  c  e 


********************************** 

ANMATJJ 


STi  (3) 
I   lT(05 
C  AprMO} 
CSiiPCft) 
CRPTLCO) 
t     SETTING    UP    TMFF 
LOl  (11 
CSwl  (1 ) 
FCl  XTf  05 
CHFCK    TO 
PF     PTCKFO 
CApO( 1  } 
OF     ANM/jTf  TP.  •> 


SAVF     THF    RETURN    ADDRESS 


MLS: 

X 

MI.SAA: 

% 
I 
% 


SAVl  J  % 

=o,?*o; 

SP3I 

SOU 

24) 
LOOP  FOR  ROW  iMpFV  TJ 

.ADRESWX    ADPFS1=A0DRFSS    OF    A*'MAT(TP.5 

6) 

SOo#M|  S2» 
SEF  WMTCH  rOS  ft  SIN  MAS  TO 

IF  JCOrl  IS  NOT  EQUAL  TO  7FR0 

501  ;  ?       TMF  criNTFNT  VF     «rl=APDRESS 

NPT  DTSTROYFD  V'lTHTN  THIS  LOOP*. 


IS 


IT 
OF 


TS  CMANGFD 
AMMATf TR.  ) 


A|  RIGHT  TO  FlNn  TMF  PPPPFR  ELEMFNTS 


X  GO  TWO 
MLS2I 

MLS1  |  t 


t    OR  COS 


I IT(2) 

ROWS  OPWN 

CAnnc 1  i 

LOAOd  -> 

STi  (?) 

cAnrM  "» 

I  0aD( 1 ^ 

STi  (?) 

LOl  (?) 

CAnD(2i 

LDI  (  3> 

LOFFl 

LOA 

LOl  (3) 

M  L  P  N 

OF  PFNPTNG 

STA 

CAnO(?'» 

LOA 

LDI  (3) 

MLRN 

ADRN 

LOl  (?) 


=  128* 

IN  ANMATrTP, ) 

FIND  TME  NEw  ADDRESS 
SAV9#0  AMD  *C0  =  0 
SAV2=CPSM  OP  COS 


SAV3rSTN  OP  SINM 
APPFS=AOnRFSS  PF  BASF 
SC2rADP.  OF  ROW  T  OF  BAE 


tc?;% 
%C?lt 
.SAV2 t% 

son 

*C?l 

.SAV3|X 

. APRFSJX 

KOJX 

.SPFCf 

SC3J 

0(?W 

.SAV?J 

SC3»X 

ON  CONTENT 

TFMP? 

sdi;x 

0(2)1 
,SAV3> 
SC3J* 
TFMPI 

.ADPFS?**  A0PFS2=AOR.  OF  TPASF 


FOR  ANMAT 


POW 
OF 


T  OF  RA^E 
.SAV2 


TtmES  COSH 


SC2=A0P.  OF  ROW  T*  1  PF  BASF 


ROW  1+1  TIMfS  SINH  OP  SIN 


00061P00 
00061 300 
000<S1  400 
00061500 
00061600 
00061700 
00061 800 
00061900 
00062000 
000621 00 
00062200 
******0006?300 

0006P400 
00062500 
0006?600 
00062700 
00062800 
0006P900 
00063000 
00063100 
O0063PO0 
00063300 
00063400 
00063500 
00063600 
O0063700 
00063800 
00063900 
00064000 
000*4100 
00064P00 
00064300 
00064400 
00064500 
00064600 
00064700 
00064800 
00064900 
00065000 
00065100 
00065200 
00065300 
00065400 
00065500 
00065600 
00065700 
00065800 
00065900 
00066000 
00066100 
O0O66PO0 
00066300 
00066400 
00066500 
00066600 
00066700 
00066800 
00066900 
00067000 
00067100 
00067?00 


-105- 


*   wr 

Ml  S3  I 


*    PF     THF 


*     .An 


TF 
SA  v 

CHA 
OF 
IF 
UP 


CAoP(?i 
STA 

I  Ol  (  ?1 
FOI  XT(?) 
I  Dl  (?) 
CCP(?1 
STl  (21 
JUST  CpANOfP 
I  PI  (?) 

papp(?i 
lpa 

S  F  r  p  N  n 
I  D I  (  3  1 

mlrm 

STa 
CAnP(?i 

LPA 
LDI  f  3  *» 

Ml  PV 
APpn 
L ni  (?i 

RFS2     IS    TMF 

CArPf?) 

f app(?} 

sta 
i  ni_  f  ?"> 

.SAV9t=0    Wf 

F     FPP    TMF     TPTAl 

N  6  E     THT     A  0  P  P  F  S  S 


scoi 

0(  ?1I 

.SAV$I 
S00#M|  S3l 
.SAV3I 

o; 

.SAV3I 
THF     SIGN    PF    SIN 
•ADRFSJ 

SCOJ 
0(?1J* 
PW    PF     F»ASF 

.SAVSJ?  *C3 
SC3J7-  POW 
TFMPI 

snu 

0(2U 
,SAV2| r 
SC  3  * 
TFMPJ 

.  appfs?; 

aodrfss  PF    TMf 
SCOJ 

spi ;*  TIN 
Of  ?  w 

.SAV^:*  CHF 

flPF    WPRtfyNr,    WyT 

N'ATPTX     ANMA 
IM     *fl     TP     P 


WF    NOW    CALCULATING    THF    Fl  FMFNTS 


=STNH    PR    -SIN. 
I     TIMFS     »C  3  I 


*C3=CPSH    PR    COS 


ANMAT(tp.1  ANP  JUMP  TO  MLSl 
.SA\'9:  =  1  WF  ARF  WORDING  WIT 
PIFFfPFNT     Fl  FMFNTS     FRpM     AnM 


TfMP.SPR.MATry  trasf 
P  RPW  1*11 

CK  ,SAV9|=0  PP  t=l. 

H  CPSM  K  StNH  WHtCH  IS  THf 

T(TR.  V  SP  WF  NFFO  NPT 

TfK    [\P    A     DIFTFRFVT    FLFMFMT 

• 

H  CPS  RSTN  AMP  HA\/F  TP  PTCK 
ATCTR.)'  Sn  WF  S"!P  TP  MLS4a 


MLS6! 
%  PFA 
MLS7: 


MLS5l 


FOl  XF(  ?1 
TXFFK«(01 

rktp 

T  X  F  F  M  f  0  5 
STAPLISH  PaSF 
I JT(01 

r apP(0> 
CPPTLf 01 

i  di  t  n 

L  p  i  ( ?  1 

PArP(  1  -» 

LOA 

r  A  p  n  ( ?  i 

STA 

TXFFMf 01 
LOl  (31 
FXTHLf 31 


SPO#MI  S6j 
#MLSl » 
*MLS7J 

*  m  l  s  n  a  r 
py  l  pap  inc. 

=0* 1 *0) 

SP31 

?«J 

.APPFS?J 

.APRFSI 

SCOJ 

0(11) 

SCOI 

0(?U 

>"LS5| 

.SAVl J 

STCRI 


TPASF  INTO  BASF 


%tt*t 

TRASP 

* 


tt1lt1t11ltttt1tttt1)l11!tt1tt*11t1tf1*t11tttttt11ttt11 

PS»     F  Tl  I   I 


STl  (31 
LOl  (01 
I  OFF! 
LOA 

i  ns 

ADM 


prpcfpupf  tra 
pf  thf  anp,!  f* 
sign  pf  thf  p 

.SAVl Xf 
.  ANTIJWJ 
SCOl 

NUMp;? 


SPPS  FINPS  TMF  TPANSPPSF 
MATPIY  BY  CHANGING  THF 
FF  -PT  AGClMAL  Fl  FMFNTS 
SAVF  PFTlJRN  ADRFSS 


FINP  THF  INPFX  FpR  THF  FIFMFNTS 
THAT  NFFP  THF  SToN  CHANGF 


*1) 
SSI 


00067300 
00067400 
O0067500 
00067600 
000^7700 
00067600 
00067900 
00066000 
00066100 
0006RP00 
00066300 
00068400 
00066500 
00066600 
00066700 
00066600 
00066900 
00069000 
000691 00 
00069P00 
00069300 
00069400 
00069500 
00069600 
O0069700 
00069800 
00069900 
00070000 
00070100 
00070700 
00070300 
00070400 
00070500 
00070600 
00070700 
00070800 
00070900 
00071000 
00071 100 
00071200 
00071 300 
00071 000 
00071500 
00071600 
00071 700 
00071800 
00071900 
O007P000 
0007?100 
0007??00 
O007?300 
00072400 
0007P500 
0007P600 
0007P700 
A007P600 
00072900 
00073000 
00073100 
00073200 
00073300 
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STA 

LPA 

CSHP(  05 

RTL 

IDFF1 

I  PA 

STA 

LDLrP5 

Lpffi 

LDS 

L  DL  (35 

I  PA 

CHSAl 

«TA       # 

LDLr?5 

FXPHl  (  3> 

*  ************* 

ANGl I     FIl I  J 

?      PRPCFM'PF  A 

*  TRAKS FORMATION 

*  SlNF=SPRT(0 
S      COST  N  F  =  S PPT 

?!     HsscPTrflYfr 

*  THrPF  APF  PITFF 
r  THF  CAI.CL'I  ATTPK' 

*  CPMF  VP     I  ki  TMF 

rtfitttttmumwit 

STi  (3) 
I  I  t  f  0  5 
CAnP(0i 
CPnT|_ro5 
LPl  f  15 
I  IK35 
L  Dl  (25 
LDrFl 

IDA 

AGt  STA 

CAPP(1 5 
TXFFM( 05 


THF  INDFX  IS  FOUND 


TFMPI 
NUMRf 

II 

SA*1  I 

scoi 

SRI 

TFMPI* 

•SPEC) 

SCOI 

TFMPI 
.ADRCI 
#0(351 


0(351 

SAV1IX     PrTuRN 

STCRI 

•♦•a*?******************************************* 


*  *  +  *  * 


mgl  ftnps  cpsinf  ano  sine*  tp  form  THF 

MATRIX  AMMAT. 

. 5-c r ?k-i  *2k1/h). 

(0.5+C C2K-1 ,?K3/H). 

?K-j,?Ki*2+(rr?K-i*?K,-n-cr?K#?k]3*2'. 

RFMT  CASFS^TP  PF  CPNSyDFRFD  FOR 

PF  CTISTME*.  THFY  WTH  PF  PlSrUSSFD  AS  THFY 
PROGRAM*. 

itmtftti*itittititiitiitittt*ittf*iiuftiit 
.SAVll*    SAVF  RETURN  APRFSS 
=0*1*01 
SD3I 
241 

.AORCI 
=  0.1 
.  S  P  F  C  I 
SC2I 

SCSI 
0(151 
SP1  I 
*  A  G  |  * 


MAKF     SI'RF    THAT    THE    PFF-DTAGO 
NAI     ELEMENTS    ARF    7FRp'. 


AG1  I 
* 


I  IT(05 
CAPP(0} 
C  S  M  P  (  0  > 

rPnTL(05 
r  L  r  ( l  S  t 

CCP(1 5 


=0*2*01 
SP3I 
SD1  I 
24j 

0(051* 


TPFATF    THF    TURN    PN    PATTFRN 
FPP    THF    OFF-HIAG.     PICC-UP 


TXFfVf 05 
STl  (15 
LDFFl 
LPA 

LlT(25 
tDv 
STA 
LDA 
RTl 

CSHP(1 5 
LDFF1 
LPA 
STA 
*  FORM  2K*2K-1  t?K' 


*  AG 
,  AN 
Sri 
NUM 
=  11 
SC2 
TEM 
MUM 
SA* 
II 
SCI 
SRI 
TFM 
1*2K 


II 
TUNI 

I 
PI 

I 

PI 
PI 
11 

I 

PI?       t  TMPEx     IS    FPU*'D    ANip 
.     T.E.  1*0*3*2*5*4*.  .,.,'. 


IT  IS  PF  THF 


00073400 
00073500 
00073600 
O0073700 
00073800 
00073900 
00074000 
000741 00 
00074200 
00074300 
00074400 
00074500 
00074600 
00074700 
00074800 
O0074900 
♦00075000 
00075100 
00075200 
00075300 
00075400 
00075500 
O0075600 
00075700 
00075800 
00075900 
00076000 
00076100 
00076200 
00076300 
00076400 
00076500 
00076600 
00076700 
00076800 
00076900 
00077000 
00077100 
00077200 
00077300 
00077400 
00077500 
00077600 
00077700 
00077800 
00077900 
00078000 
00078100 
00078200 
00078300 
00078400 
00078500 
00078600 
00078700 
00078800 
00078900 
00079000 
00079100 
00079200 
00079300 
00079400 
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*  npw  we  r  ai  r 
r  a)   HSopeov 

*  B5   H     eS 

*  C)   SlMFrcP 

*  TO  A)l 

LD« 
I.DI  ( 
I  Da 
STA 


%   appvf 


*  TO  P) I 


*  TO  C5« 


LDl  ( 
STI  f 

CLr( 
SI  tT 

cswi 

STA 

I  PI  f 

LDrF 

LDa 

MLpN 

I  iTf 
LD« 
MLRN 
STA 
LDy 
I.  DA 
I  C|  ( 
LnF 
STA 

ci  r( 

SLTT 
I  Dc 
LDTF 
LDA 
SBpN 
STA 
RTl 

rsuR 

LPrF 
LDa 
STa 
VAUtF 
LDl  ( 
LDFF 
LDA 
fcLRN 
ADRN 

C  AI  L 

STA 

I  TT( 

TAl 

SFTF 

.IAO 

SFTF 

SFTF 

LDA 
DVrN 


HI  ATF 

rr?K-i*?Ki*?*(rr?K-i*?K-n-rt?K#?K])*2 

ORT(HSOR) 
RTf0.5-Cr?K-l#?Kl/H) 

TFMPJ*     TNDFy  TN  SS  FOR  nFF-DlAGS*. 
25      .ADPCm;>   ADRFSS  OF  MATRIy  IN  USE 

#0(25J*    Cr?K-1#?K]  ! DADFP  AMP 

TFMPlj*    SAVFp  MFPF  FDR  IATFR  CASF 
DIFFFRFNTIATION  FOR  COSlNF 
3)      .NMOJ 
3)      . SAV2? 
351 

(3^  rROUTFIFXCHI  (3)  SICRJ 
(15    M 
1       SCI  J 

TFMPll?   TFMP1  HAS  IN'  FvFPY  PF  Cr2K-1.2l<i 
TN  PAIRS  OF  TWO, 
15      .SPFCl 
1       SCI  I 

TFMPl J 

TFMPlJX   £A=(C[2K-1 ,2*J5*2 
35      sfl.l 

SC3J 

$$;%  *As«y(cr?t<-i  *2ki  5*2 

TFMP2» 

MUMPJ*    TN'Pry  FDR  THr  DlAC.Fl-FMFNTS 

*0(2)  if       CC  I»  T  3  I OADFD 
15      .ANTUNJ 
F1       SCI  t 

GFRSHGJ*  SAVF  CC2K-1*2k-1] 


35f 

(35  =POUTF)FYCHL(35  *ICRI 

SA  J 
1       SCll 

GFRSHGJ 

SSJ 

GFRSHGJ?       C[2K-1  >2K-l5-Ct?lo2K1 

SA,U 
(15  D 

1  SCI  J 

SRJ 
GFRSHGJ 
S  TV  ALL  PFS  Of  GFRSHG  IN  PAIRS  PF    TWO 
15      .SPFCJ 
1       SCI  J 

GFRSHGJ 

GFPSHC,  I*       *A  =  (Cr?K-l>?K-l5-Cr?K>2i0  5*? 
TFMP?;*  SAzHSOR 


SCPT64C  ")l 

TFMP2 i* 
05      =0.1 

senj 

I  .  AND*.  FJ 
SCOJ 
j.PR.r J 

1       F.AND.F) 

TFMPl J 
TFMP2J 


H  STORED  HFRF 


00079500 
00079600 
00079700 
00079800 
00079900 
00080000 
000801 00 
00080700 
O00P0300 
000P04OO. 
O00P05O0 
00000600 
O00P0700 
OOOPOPOO 
00080900 
000P1 000 
O00P1 100 
00081200 
00081 300 
000P1 400 
O00P15O0 
00081 600 
O0081 700 
O0081P00 
00081 900 
O008P000 
00082100 
O00P??00 
O00P2300 
00082^00 
00087500 
O008P600 
00087700 
00087800 
00087900 
00083000 
00083100 
00083200 
00083300 
00083400 
00083500 
00083600 
00083700 
00083800 
00083900 
00084000 
00084100 
O00P4700 
00084^00 
00084400 
00084500 
OOOP4C-00 
00084700 
O00R4800 
00084900 
00085000 
00085100 
00085200 
O00P5300 
00085400 
00085500 
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ST* 

LDI  f  1  J 
LDrri 
i  uro) 
LDa 

SPRN 
CAI  I 
LD< 


TFMP?J* 
.SPFCI 
SCll 
sO. 51 

scot 

TFMP?t * 
SORT64C  5t 

TFMP t* 


SA«Ct?K-l/?KJ/H. 


SA  =  STNF  =  0.'i-Ct2k'-l  /2K1/W 

SS=TNDFX  FOP  C[?K/?K-ll  *  C 
pt?K-l/?K1 


LDI  f?5  .ANTUNi; 

LDFF1  SC2t 

CHOI 

LDrFl  SClt  • 

LDI  f  3>  .ADPCt 

STa  #0(35t* 

T     NOW  HEAR  THISt  CASF  DlFFF 
1!  CASE  ONF  t  C[?K-1  »?K-1  ]  - 


STNF  INTO  ANMAT 

PFNTIATIPN  FOP  COSTNF 
cr?K/?K]  LSS  0  nR  p,TR  0 


LDa 

I !Tf05 

IAI 

J  AC, 

SFTF 

SFTCfO> 

STI  (0) 

SFTF 

SFTC(0> 

C  A  K'  P  (  n  } 

7FRT(0> 

SFTF1 

1  iVf  n 

LDa 

adrm 


GFRSHftl 
=  0.  J 

scot? 
scot* 

LAND 

Ft 

.SAV9t 

J. OR.  Ft 

Ft 

SD5t* 

/AG?I* 


SA  LSS  0 
SA  GPT  0 


f; 


CLFAK'  rilT  WMFRF  *A  TOUAL  TO  7FP0 
f'P  f  r?k-l /?K-1  J-Ct?K#?Kl  LSS  OR 
CRT  7FP0/  SO  <KIP  TO  AG? 


f.ano.f; 
=0.5; 
sci  t 

TFMP?t 


CAI  I   *opT64(  )t* 

LDI (15  ,SAV9t 

LDFF1  SClt 
CH^At* 

LDfFI  SCOt 

LDY  numr; 

LDl f 25  .ADRCt 

STA  +0(?5t* 

x  cr?K-i/?K-n  -  cr?K/?Ki  l 

X*Ct*tt***ttttttttttttt1ttt 
%    CASF    TWPl     Cr?K-1 /?K-1 1    - 

AG?t     cnvpcro)} 

LDl  f  15     .SPFCt 
CA^D(M    SCOt 
7FPT(n    /AG5t* 
LDFF1       SClt* 
*  CASF  ?A»   Cr?K-l/?K]  GRT  P 


SA=S0RT(0.5+C[?K-1/?K]/H5 


SA  LSS  0  WUFRF  I  PTTS  WPF  SFT 


COSTNF  STORFD  DtAGPNALLY  WWFRF 
SS  0  OP  GRT  0 

Cr?K/?K1  FOL  0 


TF  NPNF  *RF  7FRH  FINISHED 
PFS  TURNFD  Ow  WHFPE  CASF  TWO  s  0 


I  ITC05 

IDA 

T  Ar. 

SFTF 

SFTCC35 

7FPTC35 

SFTF1 

LlTf 05 

LDa 

LDy 


=  0.t 

TFMP1 t 
SCOt 

i.and*.f; 
f; 

/AG3t* 

f.and'.f; 
=  1.; 

sco; 

numr; 


NiHNF  ARF  C,RT  0 


000R5600 
00085700 
O00R5P00 
O0085900 
00086000 
00086100 
00086200 
000P6300 
00086400 
00086500 
00086600 
00086700 
00086800 
00086900 
00087000 
00087100 
00087200 
00087300 
00087400 
00087500 
00087600 
00087700 
00087800 
00087900 
00088000 
00088100 
00088200 
00088300 
00088400 
00088500 
00088600 
00088700 
00088800 
00088900 
00089000 
00089100 
00089200 
O0089300 
00089600 
00089500 
00089600 
00089700 
00089800 
00089900 
OOOQOOOO 
000901 00 
00090200 
00090300 
00090400 
00090500 
00090600 
00090700 
00090800 
00090900 
00091000 
00091 100 
00091 200 
00091 300 
00091 400 
00091500 
00091600 
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*  CACE 
AG  3  t 


I  Dl  (?)  .ADRCI 

STA         *0f2H* 
2R|  rr?K-1 #?K]  I  SS  0 
COMPCf 3)1 


rnsiwr  ■  1  whfpf  r[?*-i»?Ki  opt  o 


CAt»P(  3} 
STI  f3) 

tnr 

I.  I  T  ( o ) 
I  Al 

SFTF 

sfttco) 

7FPT(0-) 
SFTF1 

iiTtn 

LOa 
LDv 
L ni  (?i 

STa 


FJ 


*  CASF 

AGO  I 


sen 

.SAV9I 
SC3) 
sO.I 
SCO! 
I .AND. 
FJ 

»  A  G  0  I  X 
F.ANP,FJ 
=  0.1 
$C1  » 
NUHPf 
.  AORCJ 
*0(2)  If 
0 


NPNF  ARf  |_SS  0 


2Ct   Cr?K-1*?K]  = 
CpMPCfO)* 


AG5I 

Itttttl 

H  Y  A  N  G  t 


LDl  (3) 

r  A  w  P  f  3  i 

7FPT(3} 

LPFF1 

I  iTf  1> 

LDA 

LPy 

LPl  f2) 

STA 

I  iTfl  1 

LDa 

L0< 

STA 

I  Pi  f  31 

FVrMLf 31 


.SAV9| 
SCO* 
#AG5JS 
SC3f 
=  1  .) 
SC1  t 
NUMBJ.. 
. APPCI 
♦Of  2)  J» 

=o.; 
set  J 

TFMP» 
*0C2)f 
.SAV1 »* 
MCRJ 


rnsyMFso.*    WHFPF    Cr2K-1»?Kl    LSS    0 


K'PNF    ARF    FOt     0 


rnsiNr^i .*  wHfRf  ct2K-i>2K]«o 


Pf-TllRN    TP    TMp    OMTSlPr    WPPI  D 


n  =  A 

F=A 
*2  = 

Kl  = 

TAN 
TAN 

nl  = 

Fl  = 

rPS 

r?s 

TF 
MYA 
MYA 
AT  TH 
FATTO 
FIN 


«*##**# 


Fli  I  j 

PPOrFPIlpF    HYANG    rALrilLATFSt 
C2T-1  »2*-l  J-M2T'?K]. 
[2T-1  t2Ki-A[2l,2K-l  ] '. 

siimc n   y   F)t 

SUM(n*?+r^2) , 

H=-2*2/Kl  . 

SG=TaNH*2. 

1  ./SOPTM  .-TANSp)'. 

SOPTr  f D+1 . 5/2.  J. 

HestiPTTf  F  +  l  .5/21. 

soPTr  r  (F  +  1  .  )/2]-n. 

TANH     I  S<    0     THFN     SINHe-P    F|  SF     STNP=D. 

KGt2T-lf2T-l  1  =  HYANG[?T*?n*CnSH 

KGr2T»2T-n=MYAK'Gr2I-1  »2I  1  =  STNM. 

F  FNO  PF  THF  TOTAL  CPMPllTATIPN  THF  COMvEPGFNCF 

P  IS  FPMNP  PY 

VAl  :=  o'.SYKl  Y(l  ,-SWRTORTf  1-TANRO)^. 

tttt«tt*1t*t11t1*1ft1tt1l*tttttttttttf 


S  T  i  f  3  1 

SFTF 

SFTF1 

TLRAI 

STA 

sta 


SAVl it       SAVF 

F.PR.-F) 

F.ANP.FJ 


PFTtlRW  AORFSR 


KAP1 J* 
KAP2» 


C| FAR  KAP1  AMP  KAP2 


00091 700 
O0091 800 
O0091 900 
00092000 
000921 00 
00092200 
00092300 
00092000 
00092500 
00092600 
00092700 
O0092800 
O0092900 
O0093000 
00093100 
00093200 
00093300 
O0093A00 
00093500 
00093600 
00093700 
00093800 
00093900 
00090000 
00090100 
00090200 
OO090300 
00090000 
00090500 
O0090600 
00090700 
00090800 
00090900 
00095000 
00095100 
00095200 
00095300 
00095000 
00095500 
OO095600 
00095700 
00095800 
O0095900 
00096000 
00096100 
00096200 
00096300 
00096000 
00096500 
O0096600 
00096700 
00096800 
00096900 
00097000 
00097100 
00097200 
00097300 
00097000 
00097500 
00097600 
00097700 
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HYP  I  l 


f  AnP((M 
C  St"tP  f  r>  ■> 
CPnTLC  0  5 
I  Dl  ( 15 
L  Dl  f  25 
I  DfFI 
LDy 

rs^Pd  5 

I  DFFI 

I  TTC25 

LOA 

ADM 

LOP 

LDY 

I  Dl  (  25 

LDFFl 

I  Dl  f25 

LDA 

CH«A  J 

loi  r n 

LDrFl 
STA 

LPI  (35 
STl  (35 
CLr( 35  J 
SLTT(3}    i 
LTfFI 
SBrN 
STA 

RTi 
CSuP( 1 > 

LPrrl 

STR 

ldi  f  n 
LDrri 
I  D* 
MLRN 
STA 

i  m  (35 

STI  (35 
LD« 

ClC( 351 
SLTTM}  • 
LD«; 
I  Dy 
i  Da 
CHsA  J 
IDI  f  n 

LDrFl 
STA 

r  L  r  (  3  5  j 
SLtT(35  : 
LDFFl 
SBRM 
STA 

LDI  (35 
STl  (  35 
C  L  C  (  3  5  J 


*0,2,0j*  I  nPP  FDR  FFTCPING  THF  PRPPFR  ROWS 

OF  THf  MATPTy  I'NOFP  CP*'S  IPFP  AT  ION 
SP3) 
SPU 
24J 

•  AMTllNj*  FVFPY  PTHFR  PF  TURK'FP  PN 
•SPFCJ 

SC2J 
SCOJS      TMPFX  FPP  THE  At2T-l,2K-11 
It 

sn  t 
=  1  j 

SXJ 

SC2J 
SAI 
SRJ*   INDFX  FOR  Ar?T*2«] 
.SPFCJ 
SC2J 

.APRAJ*  APRFSS  OF  PASF. 
*0(25J 

.AMTllN  J 

SCI  J 

TFMPJ?       -Ar2l-1>2K-1 ]     <TPRFD    HFRF 

.NMQJ 

.SA\/2J 

^RPtlTFJFXCHI  (35     SICRJ 

sci : 

TFMPJ?   -Ar?I*2K]-C-A[9l-l,?K-U). 

TFMPJ 

SA»1  f 

11 

SCI  J 

TFMPJ*   TFMPt=D 

.spfc; 

SCI  i 

TFMPJ 

TFMPJ 

TFMPi;*  TFMPll=D*2 

.  P  M  F  I 

.SAV2J 

SXJ 

iROMTFJFyCHI  (3)  SICRJ 
SA) 

SSJX      IVPFX  FPR  AT2I-1  »2K-T,Ar2I.2K-li 
*0(25J*  Ar?I-l»?Kl,At2T*2K-n(PADFP 

.ANTllNJ 
1  I 

SCI  J 
TFMP2J 

RntlTFIFyCHI  f3)  STCPJ 

SCI  J 
TFMP2J*   -Ar2I*2K-ll-(-Ar2T-l»2K]) 

TFMP2I 

•  NMOJ 
.SAV2I 


00097800 
00097900 
00098000 
00098100 
00098200 
00098300 
00098<iOO 
00098500 
00098600 
00098700 
00098800 
00098900 
00099000 
00099100 
00099200 
00099300 
00099400 
00099500 
00099600 
00099700 
00099800 
00099900 
00100000 
001001 00 
O0100200 
00100300 
001 00400 
00100500 
00100600 
O0100700 
0010C800 
00100900 
OOlOlOOO 
00101 100 
00101200 
00101 300 

ooiomoo 

00101500 
00101 600 
00101700 
00101800 
00101900 
00102000 
00102100 
00102200 
O0102300 
00102400 
O0102500 
00102600 
00102700 
00102800 
00102900 
00103000 
00103100 
H0103200 
00103300 
00103400 
00103500 
00103600 
00103700 
00103800 
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HYP2i  t 


HYP1  i 


SI  TTCM 

CSwl  M  t 

LPfFI 

STA 

LDl  (1  ) 

I  DffI 

IDe 

Ml  pN 

STfl 

LDa 

N  L  P  N 

ADRM 

STa 

LP* 

APRN' 

ADpM 

STfi 

TXFTM(O) 

JUMP 

I   iTfO) 
SFTF 

SFTFl 
CI  RAJ 

i  n i  tn 

l  DFFI 

IDA 

I IT(?) 

SFTF 

SFTF1 

LPc 

PTl 

t  PS 

A  PPM 

I  ITf3> 

Foi  xf  r ?■> 

ST. A 

Al  TTf  0) 
r-  P  t  P  T  f  0  ) 

JUMP 

LDl  M  ■» 
I  DFFI 

I  iWcn 

LDA 

Ml  .  PN 

PVRN 

STA 

SAP! 

LlTCO) 
TF  SO  GP  PUT  «N 
FORMATION  FOUND 

SPRM 

SAP* 

I   TTfC^ 

TAI 

SETF 

SFTTfO-) 

COMPCf 0} 

C AmP(0> 

sti  ro) 

7FRF(01 


=  POUTFIFyCMI  (3)     SICRJ 
1J 

SCI  J 

TFMP9I*    TFMP?|sF. 
.SPFCJ 
*C1J 
TEMP?) 

tf^p?) 

GFPSHGj*   r,FPSHf,jcF*2. 
TEMPJ 

tfmp?;?   n  x  f. 

KAP?  I 

KAF?J?   KAP?«=K?. 

TFMP1 » 

GFPSHGI*  P*?+F*?. 

KAP1  \ 


KAP1 J* 

hyp: 
=  0J 
E.PR.-FJ 

F. ANP.FJ 


K API  tr«1  . 


.ANT 
tfl  j 
KAP? 
si  j 

r.PR 

F.  AN 

SA»* 

$S#0 

SPj 

SS) 

*C2> 

r64J 

SC3, 

KAP? 

=  1  I 

*P1» 

HYP? 

.SPF 

sou 
=-?. 
scot 

KAP? 
KAP1 
KAP? 


II M  I 

(05) 

.'Ft 

n.ri 


LOG-SUM  FPR  KAP?  ANP  KAP1 


HYPll*   END  SUMING 

to; 

CI 


if 


KAP?trTANH 


=1.;*   CHFCK  TF  APSfTANH-n  LSS  l.P-1?. 
P  TRY  A  MFw^'ATPiy  R«SE,  FOPMEP  PY  THF  TRaNS' 
UNDER  ANGI  '. 
SCO* 

=  0.000000000001  I 
SCOJ 

FJ 


LAMP, 
El 

JD5J 
.MAXI 
*  +  ll 


oo: 
oo: 

00' 

oo: 

00' 
00' 
00' 

oo: 

oo: 

00' 
00' 
00' 

oo: 

no: 
oo: 

10' 
00' 
00' 
00 

00' 

no: 

oo: 

00' 
00' 

oo: 

oo: 

oo: 

oo: 

no: 

oo: 

ro: 

oo: 

oo: 

oo: 

oo 

oo: 

oo: 

00' 

oo: 
oo: 
oo: 
oo: 
oo: 
oo: 
oo: 
oo: 
oo: 
oo: 
no: 
oo: 
oo: 
oo: 
oo: 
no: 
oo: 
oo: 
oo: 
oo: 
oo: 
oo: 
oo: 
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JUMP 

i  nrri 
l  Da 
Ml  PN 
sTa 
CHOI 

I  iTfO) 

SP-PN 


HYP3A1 

sen 

KAPPJ 

KAP?I*   TA^?0»cTANM*? 

GFRsHO,  jScTPPF  tanso  FPR 


CAI  C.     PF    FTNVAI 


=  -1  .J 

SCO»* 


CAI I      SOPTA4C ); 


LP< 

CCPfO) 

I  Dft 

DVPN 

APRM 

I  ITCO) 

DVPM 


$A> 

sco> 

SCO! 
=  ?.J 

SCOl 


CAI I      SpPT64( )J 


JA  tr-TAM?0-C1  t  ) 
SA=SGPT(1 ,-TANSO) 

Sa :=1 ./S0PT(1 t-TANSO) 
Fl  FPUND 


I  ITfO) 
ADPN 
L  IT(2} 
PVRM 


=  1  .) 

sco> 

JC2I 


OPSH  STPPFD  TN  AnmAT 
SA  J  =  (CPSH)*? 


MNH 


*PPT64( ) I 
NUMRJ 
.  ADRC Xt 
*0(?)  J 

sa;s 
«co  j 

<  C  P  T  6  4 () J 
SAJ 
KAP?J 
=  0.1 
SCO! 

I  .and.fj 
f  .  a  n  n .  f  ; 

SSI 


SAI 
SCI  J 

kap?;?   <;tnH  fpund 

, AN TUN  I 

SC01 

NUMRJ 

=  1  J 

SC3J 

SAJ 
NUMPI 

$A»1  \ 

\\ 

SCO» 

SRJ 

SCI  I 

KAP?J 

.ADRCl 

*0C?)  J 
APF  STPRrp  PPPPFRIY.  FTNO  FtNVAI 

=  1.) 

SCO! 
GFPSMG;r    SAIsI.-GFRSHGJGFPRHGIeTANSO 


CAI 
I  Dv 
LDl  (?) 

STA 
MLPK' 

rBpn 

CAI  I 

LP* 

I  PA 

I  ITC05 

IAl 

SFTT 

SFTF1 

LDa 

CKKA  I 

LD* 

LDrFl 

STc; 

LDl  (O 
LDrFl 
LDA 

LlTf 3) 
ADM 
LD* 
LDA 
PTi 

CSPPf 0^ 
LDFF1 
L0«; 
LDrFl 
I  Da 

LDl  f?) 
STA 

ANn    c.n$v 

LlT(O) 

LDa 

SPRN' 

CAl I     S0PT64C )J 


001 

coi 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

001 

ooi 

001 

001 

001 

ooi 

001 

001 

ooi 

001 

001 

ooi 

001 

001 

001 

001 

001 

001 

ooi 

001 

001 

001 

001 

001 

001 

001 
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CH«A) 

ccpfo)  0)%  *coi«-i. 

SPRN  SCO) 

MLRM  KAP1) 

LlTfO)  =0.51 

MIrN  *COI 

STA  KAP2) 

c  l  r  f  3  )  i 

SLTTC3)  =KAP2l 

I 0AP(3>  SCOI 

STl  TO)  .FTNVAL* 

LDI  f  31  .SAVl  I 

FVrMLf?)  SICRJ 
1ttt*Htittt1l*1t1l1tttt11t*ttt1***tttf$t1l1tttf1tttf1ttttt*tt1t 
MK I     F  I  L I   t 


HYP3A  J  1 

**c# 

NULC 

%  APSCCC2K-1 ,?k 

*  If    SP>    It    tmf 

*  nf   cto.oi    amp 
y  TmF    PPOCrPURF 

*  K  T  T  H    R  P  iv     T  .     T 
It  MATRIV     PaSF*. 
tttttt*tt1tit**t*f 

STl  f  31 

i  di  rco 

LDFFl 
LDy 

I  PI  m 
I  PA 
CHsA) 

C  L  r  ( ? )  t 

I  TTf3) 


fflnnt 3> 
CSUPC31 
CRpTLC.3 
CCpf  2"> 
TXFFMC* 
ST  I  r?> 
LPrrl 
STA 

LDl  f  3^ 
STl  f  31 
CLPf 3it 
SlTTf 3> 
S  P  R  N 
I  DrFl 
STA 

LDl  f  3 ") 
STl  f31 
C.LC(  37  i 
SlTTOJ 
rSHP(?> 
I  Prrl 
STa 
LDFFl 
I  Da 
SAP!  r 
LPl  (15 
lOrFl 
TAr, 


THIS  PRPCFPURF  CHFCKSt 

-1  ]'f[?c,?Kn  t  Fo  1  .P-12', 

N  FINDS  THE  Slf,M  PF  CHM]  PppOSlTE  Pf  TMAT 
THF  RHWTK'DFV  I  =  .K'Ay.  IT  THFN  LFAPS  INTO 
SHFT  WHTPH  FyfHArQFS  THF  SEfPNp  POw 

HF  PRPCFPURF  SHFT  IS  AL SO  USFp  FPP  THF 

t*f11ttt*tt*tl9*1t1Ufft*1t*ttttt1*t*tttift 


NUl  t 


.SAVl t* 
•SPFCl 
SCO) 
NUMP) 
.  A  P  R  C  M  I 
*0( 1  )  \% 


=0>2#0** 

SD3) 
SD1) 
24) 
0(31) 
/NUl  ) 
. ANTUN) 
SC2) 
TEMP)* 

.nmp; 

.SAV2) 


SAVF  PFTtiPM  APPFSS 


LPAP  PI  AGONAL  F|  EMFNTS 


CRFaTF  PaTTfRN  TO  TURN  PN 
THF  FVFN  PFS. 


-Cf2K-l»?K-l 1  SAVFP 


:RPUTF)FyCML f 3)     STCP» 

TfMPJ*  PIFFFRFNfF    OF    OFF     DTAGS     TN    r\'FN 

SC2r*  PFS 

TEMP)?  -Ct?K#?K1-(-CC2K-l»?k-l  ]) 

.PNF) 

.SAV2) 

>RPUTF)FXCHl (3)    SICRI 
1) 
*C2l 

TFMP: 

sro» 

tfmp: 


.TMM 
SCO) 
SCI) 


) 

TW) 


TAKF  APS.  VAL.  WWERF  NFCFSS. 


oo: 
oo: 
oo 

00' 

oo: 
oo: 
oo: 

00' 

oo: 
oo: 

no: 

00' 

no' 

00' 

no 

oo: 

no: 
00 

no: 
oo: 
no: 
oo: 
oo: 

00! 

no: 
oo: 
no: 
no: 
oo: 
no: 
no: 
no' 
oo : 
no: 
no: 
no' 
no' 
oo: 
oo: 
oo: 
oo: 
no: 
no: 
no: 
oo: 
oo: 
oo: 
oo: 
no: 
oo: 
no: 
no: 

00' 

no' 
oo: 
no: 
no: 
oo: 
oo: 
oo: 

00' 
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NHL?! 


SETF 

STI  (3) 
7FpF(3l 

JUMP 

LDrFl 

LPy 

LDI  (2) 

STI   (?) 

LDA 

TSw» 

SFTF 

SF  TC(1  i 

LDI   (21 

7FrF(?1 

LEapOM W 

sktp 

I  FflHZf 1  )  J 
I IT(2) 
CAkD(?1 


STI  (2> 
FGl  XF  r?> 
JUMP 
SKTP 
ShFT!  Fill! 

STI  (3"> 
LDI  (01 
LDrrl 
I  Dl  (2) 
CAnD(2i 
LDa 

STA 
LDI  (3"> 

r  Anrc  3i 
CSl*iP(3l 

LDA 
STA 
LDA 
STA 
t  TDI  UMKS  serw  MAT 
I  iTfO) 

r AnPf Oi 
CRHTL ( 0) 
NUL3!      L  Dl  f  n 

rAnP( t  i 

LDA 

sti  (oi 
ri  rf3M 

S  L  T  T  (  3 1  = 

STA 
TVFFMf p? 

f  NOW  FXCHAMAp  srC 
LDA 

ldi  r  n 

STI  (1) 

c  l  r  r  3 )) 

SLTTf  31 


-I.anp.ei 

FJ 

. INOFXJ 
#  +  1  J? 

NUL5J 

*coi* 

NUMPF 
• ADRTMJ 
.ADPFS) 
*0(2)j* 

T . AND.FJ 
FJ 

.ENRJ 
*C1  1 
#NUI  ?»? 


NO  EXCHANGE  NFCFSSAPV  IF  TEST  TS 

SATISFIFP. 

FIND  THE  SIGN  Or  CfO*0] 


LDAD  PIAGS.  OF  CmAT 


IF  TFST  IS  SATISFIED*  CCO*01 
TS  NEGATTVF,  FLSF  PPS. 


>  +  l  J 

=  77 tfi  * 

$C\}%  TMF  TNPFX  nF  CTI.Tl  PIFFEPFNT 

IN  STGN  FPPM  P(0,0]  Is  FOUND 
.MAX) 
SD1 »+1 I 
NUL5J 
>+2J 

.  S  A  V  1  I 

tSPFf » 

SCO! 

.APRFSJ*   FXCHANGF  ROw.l  WTTH  ROW  I 

*D1  J 

0(2i; 

TEMP) 
SC2) 
SD7J 
SD1  I 
0(31  J 
0(2)) 
TFMPI 

0(3U*    FNP  FXGHAMGE  ROWS*.  Tn  EXCHANGF 
R  I  X  F  T  P  S  T  »  . 
=0.1 pOI 
*D3l 
24j 

•ADRFSJ 

sroj 

0(1  w 

.SAV2> 

routf;fxcmi  (3)  stcrj 

0(1  w 

/NUI.3I 
OND  COl.nMM  WITH  COLUMN  I 

NUMRJ 

•  ONFI 

.SAV?J 


=  ROUTF>FXCMI  (3)  tlfR) 


00122200 
00122300 
00122400 
00122500 
O0122600 
00122700 
00122800 
00122900 
O0123CO0 
00123100 
00123200 
00123300 
00123400 
00123500 
00123600 
00123700 
00123B00 
00123900 
00124000 
00124100 
00124200 
00124300 
00124400 
00124500 
00124600 
00124700 
00124800 
00124900 
00125000 
00125100 
00125200 
00125300 
001 25400 
00125500 
00125600 
00125700 
00125800 
00125900 
O0126000 
O0126100 
00126200 
00126300 
00126400 
00126500 
00126600 
00126700 
00126800 
00126900 
001 27000 
00127100 
00127200 
00127300 
00127400 
00127500 
00127600 
00127700 
00127800 
00127900 
00128000 
00128100 
00128200 
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*  rm  umns 


nui  «« 
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